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Abstract 

The complete NNLO QCD corrections to the total cross section a{e^e~ Z*,j* — > tt), 
in the kinematic region close to the top-antitop threshold, are calculated by solving the 
corresponding Schrodinger equations exactly in momentum space in a consistent momentum 
cutoff regularization scheme. The corrections coming from the same NNLO QCD effects to 
the top quark three- momentum distribution da/d\kt \ are determined. We discuss the origin 
of the large NNLO corrections to the peak position and the normalization of the total cross 
section observed in previous works and propose a new top mass definition, the IS mass Mis, 
which stabilizes the peak in the total cross section. If the influence of beamstrahlung and 
initial state radiation on the mass determination is small, a theoretical uncertainty on the IS' 
top mass measurement of 200 MeV from the total cross section at the linear collider seems 
possible. We discuss how well the IS" mass can be related to the MS mass. We propose a 
consistent way to implement the top quark width at NNLO by including electroweak effects 
into the NRQCD matching coefficients, which can then become complex. 
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1 Introduction 



It will be one of the primary goals of a future e^e~ linear collider (LC) or /i"*"^" pair collider (FMC) 
to measure and determine the properties of the top quark, whose existence has been confirmed at the 
Tevatron {Mt = 173.8 it 5 GeV Q). Although the top quark will also be object of intense studies 
at the Run II at the Tevatron and at the LHC, the measurements at a LC are important to fill the 
gaps left by the measurements in the environment of the hadron colliders. One of the most dramatic 
improvements attainable at a LC can be expected in the determination of the top quark mass. At 
the LHC, where the mass is extracted from the peak in the top invariant mass spectrum of the W 
and b originating from the top decay, a final (systematics-dominated) top mass uncertainty at the 
level of 2-3 GeV seems realistic. More precision will be difficult owing to unavoidable conceptual 
and practical problems and ambiguities in disentangling the top quark invariant mass from numerous 
effects in the environment of hadron colliders. At a LC the top quark mass can be determined from 
a measurement of the line-shape of the total cross section a{e~^e^ Z*,^* tt) for centre-of- 
mass energies around the threshold, ~ 350 GeV. The rise of the cross section with increasing 
centre-of-mass energy is directly correlated to the mass of the top quark. Because the total cross 
section describes the rate of colour singlet top-antitop events, it is theoretically and practically much 
better under control than the top quark invariant mass distribution. Because of the large top width 
(r(t bW) = ~ 1.5 GeV) the top-antitop pair cannot hadronize into toponium resonances, 

and the cross section represents a smooth line-shape showing only a moderate peak- like enhancement, 
which is the broad remnant of the IS" resonance. At the same time the top width effectively serves 
as an infrared cutoff Q and as a natural smearing mechanism [^, which allows us to calculate the 
cross section in the threshold region to high precision using perturbative QCD. It is therefore possible 
to reliably relate the cross section line-shape to the parameters of the Standard Model, most notably 
the top quark mass and the strong coupling. LC simulation studies have demonstrated that for 
50-100 fb^^ total integrated luminosity an experimental uncertainty of order 100-200 MeV can be 
expected in the top mass determination from a line-shape scan of the total cross section Evidently, 
at this level of precision an adequate control over theoretical uncertainties has to be achieved. In 
particular, a precise definition of the top mass extracted from the experiment has to be given. 

For centre-of-mass energies close to the top-antitop threshold, the top quarks are produced 
with non-relativistic velocities v <^ 1. Therefore the relevant physical scales, which govern the top- 
antitop dynamics, the top mass Mt, the relative momentum Mtv and the top kinetic energy Mfv'^, are 
widely separated. Because ratios of the three scales arise, the cross section close to threshold cannot 
be calculated using the standard multi-loop expansion in the strong coupling as, but rather a double 
expansion in and v. In the non-relativistic limit the most prominent indication of this feature 
is known as the "Coulomb singularity", which originates from the ratio Mt/{Mtv). The Coulomb 
singularity is visible as a singular (a^/f)'^ behaviour in the n-loop QCD correction to the amplitude 
7 — > it for V ^ 0. The most economic and systematic way to tackle this problem is to employ 
the concept of effective theories by using the hierarchy Mt ^ Mtv » Mtv'^ > Tt ^ ^qcd and by 
successively integrating out higher momentum effects. At leading order (LO) and next-to-leading 
order (NLO) in the non-relativistic expansion]^, the use of effective field theoretical methods seems not 
to be vital, because, clearly, in the first approximation the top-antitop pair can be described by a non- 

^ We will define what is meant by LO, NLO, NNLO, etc., in the framework of the non-relativistic expansion at the 
beginning of Sec. bl 
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relativistic Schrodinger equation [§, |5| and because, luckily, the relevant current operators do not have 
any anomalous dimension at the one-loop level. Beyond NLO, however, anomalous dimensions arise 
when relativistic effects suppressed by are included. This makes the use of effective field theoretical 
methods mandatory. In addition, the effective field theoretical approach allows for the development 
of a power counting scheme, which allows for a systematic identification of all effects contributing to 
a certain order of approximation. Those power counting rules in fact confirm that at next-to-next-to- 
leading order (NNLO) the top-antitop pair can be completely described by a conventional Schrodinger 
equation containing an instantaneous potential. In general, the same conclusions cannot be drawn for 
non-relativistic bottom-antibottom or charm-anticharm systems. 

A large number of theoretical studies at LO and NLO 0, |8|, ^ |lO| have been carried out 
in the past in order to study the feasibility of the threshold scan and other measurements at the 
top-antitop threshold. Recently, first NNLO QCD calculations for the total vector-current-induced 
cross section a{e^e~ ^ tt) have been performed in Refs. [0, 12, 13|. These analyses were 

based on non-relativistic quantum chromodynamics (NRQCD) |l5| and on the direct matching 



procedure [ffq, 17]. In this work we extend the calculations and also determine the NNLO QCD 



relativistic corrections to the top quark three-momentum distribution. Interconnection effects caused 
by gluon exchange among top-antitop decay and production processes are not considered in this work. 
They are known to vanish at NLO for the total cross section |18, 19, 2C], but can lead to sizeable 
corrections in the momentum distribution [21, 22]. We also include the total cross section and the 
three- momentum distribution induced by the axial- vector current. Because the latter quantities are 
suppressed by f ^ with respect to the vector-current-induced ones, we only determine them at leading 
order in the non-relativistic expansion. It turns out that the size of the axial-vector contributions 
is smaller than the theoretical uncertainties contained in the dominant vector-current-induced cross 
section. A discussion on the size of the axial-vector-current-induced contributions can also be found 
in ]^ ]. The three-momentum distributions presented in this work represent a first step towards an 
exclusive treatment of the top-antitop final state at NNLO close to threshold. We analyse in detail 
the origin of the large NNLO corrections to the peak position and the normalization in the total cross 



section already observed in Refs. ]11, |l^], and show that the instabilities in the peak position 
are a consequence of the use of the pole mass scheme. We show that the pole mass parameter is 
irrelevant to the peak position and define a new top quark mass, the 15* mass, which is more suitable 
to parametrize the total cross section. Whereas the 15 mass leads to a considerable stabilization 
of the peak position it does not affect the large corrections to the normalization. In this paper we 
also propose a NNLO generalization of the energy replacement rule "i? E + iTt" , by Fadin and 
Khoze, for the implementation of the top quark width by including electroweak corrections into the 
matching conditions of NRQCD. In general this leads to NRQCD short-distance coefficients that have 
an imaginary part. 

The program of this paper is as follows: in Sec. |2| we review the conceptual framework of the 
effective theories NRQCD and potential (P) NRQCD as far as it is relevant to the NNLO calculations 
carried out in this work. Section ^ contains a derivation of the integral equations that have to be 
solved and Sec. ^ describes our cutoff regularization scheme. In Sees. ^, |3| and ^ the top quark width is 
neglected. In Sec. ^ we discuss the effects of the top quark width from the point of view of (P)NRQCD. 
Some details about our numerical methods to solve the integral equations are given in Sec. |6[ A first 
analysis of the total cross section and the three momentum distribution in the pole mass scheme is 
given in Sec. 0. Section ^ concentrates on the origin and interpretation of the large NNLO corrections 
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in the pole mass scheme and introduces the IS* mass. The relation of the IS mass to other mass 
definitions is discussed. Section ^ contains our conclusions. In Appendix |A|, details of the NNLO 
matching calculation are given in the framework of our regularization scheme. 



2 The Conceptual Framework - NRQCD and PNRQCD 



In this section we review the effective field theories NRQCD and PNRQCD, which form the conceptual 
framework in which the NNLO corrections to the top-antitop cross section close to threshold are 
calculated. By N'^LO {k = 0, 1, 2, . . .) for the total cross section we mean a resummation of all terms 
proportional to a^f", with m + n = + in perturbation theory in as supplemented by a 

subsequent expansion in the top quark velocity, i.e. in the limit ^ f <C 1. Thus at the NNLO 
level all terms proportional to vJ2'^=o{'^s/'v)'"'[^',as,v;a'l,asV,v'^] have to be resummed to all orders 
in conventional perturbation theory in Og, where the dominant terms in the non-relativistic limit are 
determined by a non-relativistic Schrodinger equation with a Coulomb potential V{r) = —Cpas/r, 
Cp = 4/3. In this context one has to count the strong coupling ag of order v, as long as the 
renormalization scale is much larger than the typical hadronization scale Aqcd- For simplicity we 
postpone the effects of the top quark width until Section |5|. 

NRQCD jlj, |l5[ is an effective field theory of QCD specifically designed to handle non- 
relativistic heavy quark-antiquark systems. NRQCD is based on the separation of the low momentum 
scales MtV and MtV^ , which govern the non-relativistic quark-antiquark dynamics, from the high mo- 
mentum scale Mf , which is relevant for hard effects involved in the quark-antiquark production process 
and quark-antiquark and quark-gluon interactions. The NRQCD Lagrangian for the top-antitop sys- 
tem is obtained from QCD by integrating out all hard quark and gluon momenta of order Mf or larger, 
and the corresponding antiparticle poles of the small components. Treating all quarks except the top 
as massless, the resulting non-renormalizable Lagrangian reads 

^NRQCD = - ^ Tr G^'^ G^, + 'i^Pl 

q=u,d,s,c,b 



iDt + C2 



2Mt 



+ C4 



8 Mi 



+ 



CF9s „ , cogs 



2Mt 



8Mf 



D-E-E-D) + ^icT{DxE 



E X D) + 



2Mt 



C4 



8Mf 



+ 



CFQs 



a B + 



cp gs 
8 Mi 



D E - E D) + 



csgs 



ia (D X E - E X D) + 



X: 



a) 



2 Mi 8Mf ' ' 8MI 

where only those terms are displayed explicitly which are relevant to the NNLO calculations in this 
work. The gluonic and light quark degrees of freedom are described by the conventional relativistic 
Lagrangian, where G^^ is the gluon strength field tensor, q the Dirac spinor of a massless quark. 
The non-relativistic top and antitop quark are described by the Pauli spinors V' and Xi respectively. 
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For convenience all colour indices are suppressed and summations over colour indices are understood. 
Dt and D are the time and space components of the gauge covariant derivative Dfj,, and = 
and = ^e^^^G^^ the electric and magnetic components of the gluon field strength tensor. The 
short-distance coefficients C2, C4, cp, cd,cs are normalized to one at the Born level. The subscripts F, 
D and S stand for Fermi, Darwin and spin-orbit. We emphasize that the mass parameter Mt used 
for the formulation of NRQCD is the top quark pole mass. Although it is known that this choice can 
lead to a bad behaviour of the perturbative coefficients at large orders, the pole mass is still the most 
convenient mass parameter to be used at this stage, because the formulation of NRQCD is particularly 
simple in this scheme. 

In addition to integrating out hard quark and gluon momenta and the small components in 
the QCD Lagrangian one also has to do the same in the vector [j""^^ = t^^jt) and the axial-vector 
currents (j^ = ^7^75^), which produce and annihilate the top-antitop pair close to the threshold with 
centre-of-mass energy \f(p. This means that we have to expand the respective QCD currents in terms 
of NRQCD currents carrying the proper quantum numbers. In momentum space representation the 
expansion of the QCD vector current in terms of NRQCD currents reads {k = 1, 2, 3) 

= c^(v;W)(9)-g|^(^V;t(-iB)2x)(g)-f..., (2) 
~3U-Q) = c^(xW)(-g)- A(xW-iBf^)(-g) + ..., (3) 



t 

and the expansion of the QCD axial-vector current in terms of ^Pi NRQCD currents 

jUq) = ^(v^n-|BxcT),x)(g) + ..., (4) 

~fk{-<l) = §-^{x\-lDxa%^){-q) + ..., (5) 

where the constants cf 2 and cf are the short-distance coefficients normalized to one at the Born level. 
The time components of the currents do not contribute because the trace over the massless lepton 
fields that describe the e+e" annihilation process is proportional to {S^^ — e'e-'), {e^,e^,e^) being the 
unit-vector pointing into the centre-of-mass electron direction. In addition, the zero component of the 
vector current vanishes. The dominant NRQCD current in the expansion of the QCD vector current 
has dimension three. Thus for a NNLO description of the cross section we have to expand the QCD 
vector current in NRQCD currents up to dimension five. The QCD axial-vector current only needs 
to be expanded up to dimension four. For the NNLO calculation of the cross section only the 0(0^) 
short-distance corrections to the coefficient c\ have to be calculated. 

To formulate the total tt production cross sections in e+e" annihilation in the non-relativistic 
region at NNLO in NRQCD, we first define the vector and axial-vector-current-induced cross sections 
using the optical theorem and starting from their corresponding expressions in full QCD: 



-i J d^xe'i-'^ {0\T jl{x) f\0)\0) 



^ ^l^[-i{0\TjUq)f^{-q)\0)], (6) 



R'^{q^) = ijim 



-i J d^xe''^-^0\Tjf{x)riO)\0) 
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= ^lm[-i{0\Tj'^{q)~f\-q)\0)]. 
In terms of R" and i?" the total cross section cr^^f {e~^e~ — > 7*, Z* tt) reads 



(7) 



Ve Vt Qt + 



f-Ml 











2 1 2 





q^-Ml 



where 



Vf 

a/ 



2 I 2 



3g2 ' 

tI -2Qf sin^gyy 
2 sin 0w cos 9 



2 sin Qy^ cos 



(8) 

(9) 
(10) 
(11) 



Here, a is the fine structure constant, Qt = 2/3 the electric charge of the top quark, 6w the Weinberg 
angle, and refers to the third component of the weak isospin; QtK" is equal to the total normalized 
photon-induced cross section, which is usually referred to as the i?-ratio. To determine K" and R"" at 
NNLO in NRQCD we insert the expansions in Eqs. d)-® into Eqs. (|) and (0). This leads to the 
expressions 



-n-NNLOW ) 

-Rnnlo(9^) 



47r 



where 



A"" 
A" 



i( 



6M2 



Im 



A\q^) 
A^^q^) 



+ ... 



+ ... , 



1 



6M/ 



i( 



and 



1. 



(12) 
(13) 

(14) 
(15) 

(16) 
(17) 



The expressions for i?" and R°^ at NNLO in the non-relativistic expansion in Eqs. (|T2|) and (13) 



represent an application of the factorization formalism proposed by Bodwin, Braaten and Lepage |15|| . 
The cross sections are written as a sum of absorptive parts of non-relativistic current correlators, each 
of which is multiplied by a short-distance coefficient. The vector correlator A" describes the top- 
antitop system produced in an S-wave spin triplet state and the axial-vector correlator A"' describes 
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the system in a corresponding P-wave triplet state. The axial vector-current-induced cross section is 
suppressed by with respect to the vector-current-induced cross section. We note that the correlators 
are defined within a proper regularization scheme. For convenience we have not expanded into 
a sum of a dimension six and a dimension eight current correlator. The short distance coefficients 
C" and encode the effects of quark and gluon momenta of order the top mass or larger in top- 
antitop production and annihilation vertex diagrams. As the non-relativistic correlators they are 
regularization-scheme-dependent. In principle the non-relativistic correlators and A"' could be 
calculated from the Feynman rules derived from the NRQCD Lagrangian (|l]). Such a task, however, 
would be quite cumbersome, because an infinite number of diagrams would still have to be resummed. 
Clearly, we would like to derive an equation of motion for the off-shell top quark four point Green 
function, which in the non-relativistic limit reduces to the well known Schrodinger equation, which 
automatically carries out the resummation of the relevant diagrams. A formal and systematic way 
to achieve that is to also integrate out top quark and gluon modes carrying momenta of the order 
the inverse Bohr radius ~ Mtv. The resulting effective field theory is called "potential NRQCD" 
(PNRQCD) [^4|. The basic ingredient to construct PNRQCD is to identify the relevant physical 
momentum regions in the description of heavy quark-antiquark systems in the framework of NRQCD. 
Those momentum regions have been identified in Ref. p^ j by constructing an asymptotic expansion of 
Feynman diagrams describing heavy quark-antiquark production or annihilation close to threshold in 
terms of the heavy quark centre-of-mass velocity. Because NRQCD is not Lorentz-covariant, the time 
and the spatial components of the momenta have to be treated independently. There are momentum 
regions where time and spatial components are of a different order in the velocity counting. The 
relevant momentum regions are "soft" {k^ ~ Mtv, k ~ Mfv), "potential" (k^ ~ Mfv'^, k ~ Mfv) 
and "ultrasoft" (A;'' ~ Mfv'^, k ~ Mtv^). It can be shown that heavy quarks and gluons can have 
soft and potential momenta, but only gluons can have ultrasoft momenta. A momentum region with 
k^ ~ Mtv, k ~ Mtv^ does not exist. It is in principle not excluded that there are momentum 
regions scaling like M^f " with n > 2, but even if such regions existed they would be irrelevant to top 
quark production because they would represent modes below the hadronization scale. To construct 
PNRQCD one integrates out "soft" heavy quarks and gluons {k^ ~ Mtv, k ~ Mtv) and "potential" 
gluons (A;'' ~ Mtv^ , k ~ Mtv), supplemented by an expansion in momentum components of order 
Mtv^ . Heavy quarks carrying potential momenta and gluons with ultrasoft momenta are kept as 
dynamical fields. This leads to spatially non-local four (heavy) quark operators which represent a 
coupling of a quark-antiquark pair separated by distances of order of the Bohr radius ~ 1/Mtv. For 
a quark-antiquark pair in a colour singlet state this non-local interaction is nothing else than the 
instantaneous potential of a quark-antiquark separated by a distance of order the inverse Bohr radius. 
Generically the PNRQCD Lagrangian has the form 

■^PNRQCD = -^NRQCD ^ j ^ ^ 

where the tilde above >Cnhqcd on the RHS of Eq. ( p^ ) indicates that the corresponding operators 
only describe potential quark and ultrasoft gluonic degrees of freedom. In addition, an expansion in 
momentum components ~ Mtf ^ is understood. V is the heavy quark-antiquark potential and is given 
below. Using the velocity counting rules for potential heavy quarks mentioned above we see that the 
LO contribution to V , the Coulomb potential —CFOis/\r\, counts as , i.e. it is of the same order 
as the kinetic energy. Thus, as is well known, in the non-relativistic limit, the Coulombic interaction 
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between the heavy quark pair has to be treated exactly rather than perturbatively. We note that 
the expansion of the heavy quark currents in terms of NRQCD currents (Eqs. )~(|lt)) is in general 
also affected by going from NRQCD to PNRQCD, because also in the NRQCD currents soft and 
potential gluonic degrees of freedom have to be integrated out. However, at NNLO this does not 
affect the results displayed in Eqs. @ and (P). The only (and fortunate) practical consequence is 
that we can neglect the gluonic contribution in the covariant derivatives. Collecting all terms from 
the PNRQCD Lagrangian which contribute at NNLO, i.e. count as f^, or v'^, one can derive 
the following equation of motion in momentum space representation for the Green function of the 
time-independent Schrodinger equation, valid at NNLO in the non-relativistic expansion: 



Mt 



where 



4M3 





P'o \] 




AMf ) . 




V{k,k') 



G{k,k';q^) + J ^V{k,p')G{p',k';q') = {2 7rf S^'^k 



(2vr) 



Vc{k - k') + FBF(fc, k') + V^^{k - k') , 



k'), 
(19) 

(20) 



and 



pI 



(21) 



The parameter pQ is equal to the centre-of-mass three- momentum of the top quarks. We have chosen 
this rather unusual representation for the energy parameter since it greatly simplifies, because of its 
symmetric form in Eq. (^), the analytic matching calculations we carry out in App. The same 



trick has already been used in [17|] and later in 1 11, 12, 13 1. It is a non-trivial fact that the PNRQCD 
operators, which describe interactions of the heavy quarks with ultrasoft gluons, do not contribute 
at NNLO, and that a conventional two-body Schrodinger equation with an instantaneous potential is 
fully capable of resumming all terms oc a^v"' that belong to NNLO. We come back to this issue at 
the end of this section. 

The individual potentials in momentum space representation read (a^ = as(//), Ca = 3, 
Cf = 4/3, T = 1/2, Q = k- k'): 



Vc{k) 



VBF{k,k') 



At: Cp CLs 

fc2 



1 + 



+ ai 



2 r 



StSi 



(2/3oai In(^) +02 

[QSt){QSi 



(22) 



M2 



{k + k'f 



{k' - k") 



nCpas {St + St){k x A;') 



(23) 



V^NA(fc) 



TT^ Ca Cf al 
MAk\ 



(24) 
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where St and are the top and antitop spin operators and (n; = 5) 
Po = Y 3 

Pi = Y^A-fcATni-ACFTni, 

31 ^ 20 

ai = -CA--Tni, 

f 4343 o TT^ 22 \ o / 1798 56 \ ^ ^ 

y-16C3)Ci.rn;+(^^rn;)'. (25) 

The constants (3q and /?i are the one- and two-loop coefficients of the QCD beta function and 
7 = 0.577216 ... is the Euler constant; Vc is the Coulomb (static) potential. Its 0{as) and 0{a1) cor- 
rections, which come from loops carrying soft momenta, have been determined in p6| , |27[ | and p8| , |29[ |, 
respectively.0 Vbf is the Breit-Fermi potential known from positronium. It describes the Darwin and 
the spin-orbit interactions mediated by longitudinal gluons and the hyperfine interactions mediated 
by transverse gluons in the potential momentum region. Then, I^a is a purely non-Abelian poten- 
tial generated (in Coulomb gauge) through a non-analytic term in a soft momentum one-loop vertex 



correction to the Coulomb potential involving the triple gluon vertex |3C, ^ (see also Ref. 

We have already noted that it is a non-trivial fact that ultrasoft gluons do not contribute at 

NNLO, which means that a common two-body Schrodinger equation, i.e. a wave equation containing 

an instantaneous interaction potential, is indeed capable of resumming all the terms that we count 

as NNLO. Although this is a well accepted fact for positronium and the hydrogen atom in QED p^ , 

it is not at all trivial to understand, even in the QED case, if one has to rely on arguments that 

are not in the framework of (P)NRQCD. In QED the corrections caused by ultrasoft photons are 

known as retardation effects. The Lamb shift in hydrogen is the most famous example. Based on 

the identification of the various momentum regions for quarks and gluons mentioned above, however, 

one has transparent power counting rules at hand; these show that the non-instantaneous exchange of 

gluons among the top quarks does not lead to any effects at NNLO. For the validity of the argument 

it is important that the scale Mtv^ is much larger than the typical hadronization scale Aqc^).^ To 

see that retardation effects are suppressed by at least three powers of v with respect to the non- 

relativistic limit (LO), we recall that the only source of non-instantaneous interactions in PNRQCD 

are the ultrasoft gluons, for which all momentum components are of order Mtv^. In addition, only 

transversely polarized gluons need to be considered as ultrasoft, since we can work in the Coulomb 

gauge where the time component of the longitudinal gluons vanishes. Thus, an exchange of an ultrasoft 

gluon among the heavy quark-antiquark pair is already suppressed by with respect to LO from the 

coupling of transverse gluons to the heavy quarks. To see that an additional power of v arises from 

the loop integration over the ultrasoft gluon momentum, we compare the f-counting of the integration 

measure and the gluon propagator for ultrasoft and potential momenta. In the ultrasoft case, the 

^ The constant a2 was first calculated in Ref. |^^. In Ref. an error in the coefficient of the term oc vr^C^ was 
corrected. 

Even for energies very close to threshold, the scale Mtv^ cannot become smaller than Ft since the dominant effect 
of the top width is to effectively shift the energy into the positive complex plane by an amount Ft (see Section |^ . 
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product of the integration measure d^/c and the gluon propagator counts as xv~'^ = v'^, whereas 
in the potential case the result reads x v~'^ = v"^. Because the potential momenta contribute at LO 
we find that the ultrasoft gluons can indeed only lead to effects beyond NNLO. Even if the gluon self 
coupling is taken into account, this conclusion remains true, because it also leads to additional powers 
of V if ultrasoft momenta are involved. We note that in arbitrary gauge the conclusions are true only 
after all gauge cancellations have been taken into account. The relation Mtv"^ ^ ^qcd is needed for 
the above argumentation: otherwise, the coupling of ultrasoft gluons to the heavy quarks or among 
themselves, as{Mtv^), could be of order 1. Therefore our conclusion that retardation effects do not 
contribute at NNLO would not be valid for the bottom quark"^ or even the charm quark case. 



3 Lippmann-Schwinger Equation 

The non-relativistic current correlators in the NRQCD factorization formulae for the total top-antitop 
cross section close to threshold, Eqs. ( |T2[ ) and (13), are directly related to the Green function 
G{k,k';q'^) of the Schrodinger equation (|l^), which describes off-shell elastic scattering of a top- 
antitop pair with centre-of-mass three momentum ibfc into a top-antitop pair with three momentum 
zizk' . We emphasize that the Green function does not describe the scattering of on-shell top quarks 



because the common three-dimensional formulation in the form of the Schrodinger equation (19) al- 
ready contains an implicit integration over the zero-components of the momenta in the heavy quark 
propagators. Because the heavy quark potential is energy-independent this integration is trivial by 
residues. In the first part of this section we give the relations of the current correlators to the Green 
function in the three-dimensional formulation. In the second part, we present the generalization to 
four dimensions for those results that are essential for a proper treatment of the W~^W~bb phase space 
once the top quark width is taken into account. For this section we still assume that the top quark is 
stable. 

Taking into account the partial wave decomposition of the Green function of the Schrodinger 
equation ( p^ ) 

oo 

Gik,k') = Y.G\k,k'), (26) 

1=0 

where I is the total angular momentum quantum number, we find the following relation between the 
NRQCD non-relativistic current correlators ( [l4[) and (|l5| ) and the S and P wave contributions to G: 



* For the case of bottom-antibottom quark sum rules the conclusions can, however, still be correct if the effective 
smearing range is chosen larger than P^qcd- [0 
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(fk f Sk' fcxcr~i,, ,,.k'xcr 
' -G (k, k ^ 



(27r)3 J (27r)3 Mt 
Sk f Sk' kk' ~i 



(27r)3 J (27r)3 



2 ^ (^1 ^ ) 



(27) 



where a proper UV regularization is understood. For simplicity we have dropped the energy argument 
of the Green function G. In this work we use the Lippmann-Schwinger equation, the Fourier transform 
of Eq. ( [TqI ) with respect to fe', 

^ ' ' ' G{k,x) = exp(ifca;) 



(pI 


pI \] 




mf ) _ 



2 7r 



V{k,p')Gip',x) 



(28) 



to derive integral equations for G^ and G^, which are then solved numerically. Using the partial wave 
decomposition of exp(ikx) (k = \k\, x = \x\): 



exp{ikx) = exp{ikxcosd) = y^ i'' {21 + 1) ji {k x) Pi (cos 6) , 



(29) 



where ji are the spherical Bessel functions ji{x) = and Pi the Legendre polynomials, 

one arrives at the following equations for G^ and G^: 



k' 





Mt AM? 



Po 



4 

Po 



Mt 4M? 



Mt 



Po_ 

Mt 



G°ik,x) 
G\k,x) 



sm[K X ) 
k X 



d^p' 
(2^ 



V{k,p')G\p',x), (30) 



3i 



kx / s'm{kx) 



k"^ x"^ 
d^p' 



cos(/c x) 



k X 

V{k,p')G\p\x). 



(31) 



Because in Eq. ( pO|) only S-wave states are considered, one can, instead of the complicated form of 
Vbf) use the angular average with respect to the angle between p' and k of the 1/Q^ term on the 
RHS of Eq. (p3[). Evaluating also the spin matrices for the S-wave state, the Breit-Fermi potential 
simplifies to 

1 



vi^{k,k') = — J dnv^^{k,k') 

11 vr Gf as 2 ^ '^s + ^ 



/2 



(32) 



In Eq. (|3l|), on the other hand, we will just use the LO Coulomb potential because the axial- vector 
contribution to the total cross section is already suppressed by f^. For the same reason we do not 
include any kinematic relativistic corrections in Eq. (^T|) . Defining the S-wave and the P-wave vertex 
Green function as 



Sik) 
P{k) 



d'p' ^0 



(2vr)= 



G\k,p') 1 + 



6M2 



''P' ^G\k,p') 



(2 7r)3 fc2 



(33) 
(34) 
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we finally arrive at the integral equations for S{k) and P{k), 



Sik) 



Gf{k) ( 1 + 



Gf{k) 



(2^ 



V{k,p')Sip') 



P{k) = (k) - G^ (k) J 



^^yrik,p')Pipr 



where 



G^ik) 



k — — ie 



1 + 



AMI 



(35) 
(36) 

(37) 



is the free vertex function. We note that (27r)^ 5^^\k — k') G^ (k) is the Green function of Eq. (|T^) for 
V = 0. 

The vertex functions S{k) and P{k) only depend on the spatial momentum k. As mentioned 
at the beginning of this section, their dependence on the time component has been eliminated by 
a trivial integration by residues. For a proper integration over the phase space of the top-antitop 
decay products, which is carried out in Sec. ^, we also need the full dependence of the vertex functions 
on k^. It can be recovered by comparing G-^ {k) to the product of a free PNRQCD top and antitop 

2 4 ,24 

propagator at NNLO carrying momenta kt = (A:*^ + — gff^))^) ^^'^ % = (^^ ~ (jM^ ~ sMf-''^-^' 
respectively^: 

G^{k^,k) 



k^ + 



'jl \ _ ( 

-2Mt SMf > ^ 



^2Mt SM'^^ 



k' 



k^ 



2Mt 8Mf 



le 



)L02 _ (n 

\ 2Mt 2M-, 



+ 



\ 2Mt 2Mt 



2 \ 2 



7 pI , e 

\ 2Af 2 "I" 2Af2 



kL 

2Mt 



+ ^e 



2 • 



(38) 



We emphasize that an expansion of the NNLO relativistic effects is understood. The relation between 
Gf{k) and Gf{k^,k) reads 



+00 

Gf{k) = -^ j '^G^ik^k). 



(39) 



Recalling that the potentials as well as the production and annihilation vertex corrections do not 
depend on the zero-components of the momenta, the integral equations for the generalized vertex 
functions read 

5(fc°,fc) = Gf{k',k)l^l + -^^ +^Gf{k',k) J^/j ^V{k,p')S{p",p'), (40) 



P{k^,k) = Gf{k^,k) + iGf{k^,k) J-^ J 



f dp' ^pv-(fc,p')^'(p'°,p'). 



27r 



(41) 



^ This choice corresponds to a situation in top-antitop ladder diagrams where half of the centre-of-mass energy is 
flowing through the top and half through the antitop line, see Fig. |l|. 
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The three- and four-dimensional versions of the vertex functions S and P are related by an equation 
similar to Eq. (^). We note that, by construction, the relation 



X(fc°,fc) 



X{k) 
Gf{k) ' 



{X = S,P), 



(42) 



holds for the amputated vertex functions. 

It is straightforward to formulate the optical theorem, which relates the imaginary part of 
the correlators A" and A"" to explicit phase space integrals over the modulus squared of the vertex 
functions S and P. The relations read 



Im 



12 vr^ Nc 
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+ k^ 



(43) 
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.2Mt 


2Mt 



P{k^,k) 

Gf{kO,k) 



k' 



Po_ 

2Mf 2Mt 



d^fc k^ 



(27r)3 M2 



Gf{k) 



2 Jk^ 



pI 



\Mt Mt 



(44) 



4 Regular izat ion Scheme and Short-Distance Coefficients 

All equations derived previously have to be considered in the framework of a proper UV regular- 
ization scheme. In fact, UV linear and logarithmic divergences arise in Eq. ( |35| ) from the NNLO 
non-Coulombic potentials and from the kinetic energy and vertex corrections. However, we emphasize 
that even in the case when no UV divergences arise, all integrals have to be consistently regularized, 
because only in a consistent regularization scheme can the short-distance coefficients be defined prop- 
erly. In principle, the regularization scheme "of choice" would be an analytic scheme like MS as it 
is usually used in modern perturbative QCD calculations. The preference for the MS scheme arises 
from the fact that it naturally preserves gauge invariance. Ward identities and, particularly important 
in the framework of effective field theories, power counting rules. Unfortunately an analytic solution 
of Eqs. (|T9|), ( |35[ ) and ( p6|) is not even known for four dimensions. Thus, the only sensible way to 
use the MS scheme is to start from the known Coulomb solution of the non-relativistic Schrodinger 
equation and include NLO and NNLO corrections via time-independent perturbation theory, and then 
explicitly construct the spectral representation of the Green function at the NNLO level. While this 
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(q, 0) 





Figure 1: Routing convention for loop momenta in ladder diagrams. 



program might still be feasible for the determination of the total cross section, it is rather cumbersome 
for the calculations of distributions at NNLO. 

In this work we use a momentum cutoff regularization scheme by simply excluding momenta 
that have a spatial component larger than the cutoff A. This is in fact the most natural regularization 
scheme for a numerical solution of Eqs. ( [35| ) and (36). However, a cutoff scheme contains a number of 



subtleties, which shall be briefly discussed in the following. As indicated before, a cutoff scheme leads 
to violations of gauge invariance and Ward identities in the (P)NRQCD calculation. These effects, 
however, are generated at the cutoff and are therefore cancelled by corresponding terms with a different 
sign in the short-distance coefficients. Thus the cross section, which contains the proper combination 
of non-relativistic correlators and short-distance coefficients, is gauge-invariant and satisfies all Ward 
identities up to terms beyond the order at which the matching calculation has been carried out. An- 
other subtlety of a cutoff prescription is that it is only well defined if a specific routing convention for 
the momenta in loops is adopted. For the calculations of the top-antitop cross section at NNLO, it 
is straightforward and easy to find such a routing convention, because only ladder-type diagrams are 
involved. It is natural to choose the routing used in the integral equations ( |35| ) and (^), which we 
have, for clarity, depicted graphically in Fig |l|. It is important to use exactly the same routing for the 
(P)NRQCD diagrams calculated in the matching procedure to obtain consistent results for the NNLO 
short-distance coefficients. Finally, it has to be mentioned that a cutoff scheme inevitably leads to 
power counting breaking effects. This means in our case that a term in the Schrodinger equation (|l9|), 
which is NNLO according to the velocity counting, can in principle lead to lower order contributions in 
the non-relativistic current correlator. Like the terms that violate gauge invariance and Ward identi- 
ties, the power counting breaking terms are also generated at the cutoff; they are therefore cancelled in 
the combination of the correlators and the corresponding short-distance coefficients. However, all this 
happens only if the cutoff is chosen of the order of Mj . To illustrate this issue let us first consider the 
LO two-loop NRQCD vector current correlator, which contains the exchange of a Coulomb gluon as 
shown in Fig. |^a. It is straightforward to calculate the absorptive part of the diagram as an expansion 
in pq (see App. 



(1) _ Cf as Afj 



(45) 



The first term in the brackets on the RHS of Eq. ( pS] ) is the well known Coulomb singularity, which 
leads to a finite cross section at order a^. The second term is cutoff-dependent and leads to a short- 
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a 



b 



Figure 2: The two-loop NRQCD vector current correlator with the exchange of one Coulomb gluon 
without relativistic corrections (a) and with the kinetic energy corrections (b). 



distance correction oc (In App. ^ the reader can convince her /himself that the scale of Us in 

this term is indeed of order M^.) Obviously, to avoid a breakdown in the separation of long- and 
short-distance contributions^, A has to be chosen of order Mj. (Choosing A of order Mtv would be 
absurd anyway, because this would cut off a large part of the dynamics. In this respect the cutoff A 
acts in a way completely different from the "cutoff" scale in the MS scheme, which would naturally 
be chosen of order Mtv.) Let us now consider the kinetic energy corrections in the same two-loop 
NRQCD diagram, as shown in Fig. |2|b. Using the result obtained in App. ^ the expression for the 
diagram including combinatorial factors reads 



r(i) 



Mi 



47r2 



Mi 



+ 



2 Mi 



(46) 



As expected from power counting arguments, the kinetic energy correction leads to a contribution 
suppressed by p^/Mf with respect to the pure Coulomb exchange diagram in Eq. (45). However, 



there is also a term proportional to po/Mt because A is of order Mj. This is an example for a power 
counting breaking term. In the short-distance coefficient this term leads to a NNLO contribution 
oc a^j-^. Because this term only arises if NNLO relativistic effects are taken into account, we have 
to count it as NNLO. Similar terms are caused by the Breit-Fermi potential Vbp, the non-Abelian 
potential V^a and the dimension-5 NRQCD vector current. We emphasize again that, as is the case 
for the terms that violate gauge invariance and Ward identities, all power counting breaking terms are 
automatically cancelled in the combination with the short-distance coefficients up to terms beyond 
the order at which one carries out the matching procedure. In our case, where the matching is carried 
out at order a^, power counting breaking terms of order remain uncancelled, but they are beyond 
NNLO accuracy. 

Taking into account the issues discussed above, it is straightforward to determine the matching 
coefficient C^. Details of this calculation are given in App. The result reads (a^ = as(/i)) 



1 + 



4C_F as 



vr 



Mt 



NLO 
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SttM, 7r2 
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Is? + 12 Mi _^ A^ - 12 Mi 



6AMt 



6 AM. 



+ 2 ln( — 



In this context "long-distance" efltects are not understood as "non-perturbative" effects, which come from scales of 
order Aqcd , but rather as effects governed by scales of order Mtv or Mtv^ . 
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(47) 



As explained in App. ^ we have displayed the NLO and NNLO contributions to C" separately. 
Appendix ^ also contains a discussion on the convergence of the perturbative series in in our 
cutoff scheme compared to the MS scheme [35, 36 1. 



5 Top Quark Width 

Up to now we have treated the top quark as a stable particle. For top quark production close to 
threshold this is not appropriate because the kinetic energy of the top quark ~ Mtv'^ ~ Mto^ ~ 2- 
3 GeV is of the same order as its decay width. It has been shown by Fadin and Khoze ^ that 
in the non-relativistic limit the top width can be consistently implemented by calculating the cross 
section for stable top quarks supplemented by the replacement E ^ E + iVt , where E is the centre- 
of-mass energy measured with respect to the two-particle threshold. For the total cross section, this 
prescription remains valid even at NLO because, in this case, order as QCD radiative corrections in 
form of a gluon that connects top quark production and decay vanish ||2^ . A consistent implementation 
of the top width at NNLO has been missing so far. 

In the framework of NRQCD, and if one is not interested in any differential information on 
the top decay products, the top width can be understood as a modification of the NRQCD matching 
conditions caused by electroweak corrections. Because the particles involved in these corrections can 
be lighter than the top quark (i.e. if the top quark can decay weakly) they can lead to non-zero 
imaginary parts in the matching conditions and, likewise, in the short-distance coefficients of the 
NRQCD Lagrangian and the NRQCD currents. This is a well known concept in quantum mechanics 
of inelastic processes where particle decay and absorption processes are represented by potentials 
and couplings carrying complex coefficients, if one is not interested in the details of the decay and 
absorption process. In this context the effects of the top quark decay are only a small (but nevertheless 
the most important) part of a whole array of electroweak corrections relevant to top-antitop quark pair 
production close to threshold. In this work we only consider the effects from the on-shell top decay 
into a W boson and a bottom quark, assuming that the W and the b are themselves stable. Thus the 
results presented here are, by definition, gauge-invariant. The consistent treatment of all electroweak 
effects, including a proper handling of the off-shell decays of the top quark, interconnection effects 
and of gauge-invariance-violating contributions, is beyond the scope of this paper. Such a treatment 
is carried out in [37|. 



Before we determine the modifications of the NRQCD Lagrangian, Eq. (jl]), through the top 
decay, some remarks about the velocity counting of the top quark width are in order. Comparing the 
numerical size of the top width 
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with the binding energy of a fictitious toponium IS" bound state in the pole mass scheme (a^ ~ 0.15) 



1.75 GeV. 



(49) 



we find that we have to count Tt/Mt as order , i.e. 



Mt 



1 2 



(50) 



Recalling the velocity counting of the operators of the NRQCD Lagrangian, this means that at NNLO 
and to first order in Tt/Mt we have to determine the coefficients of the NRQCD operators i/j^iFttl^, 
il^Hjf^Dttl^, ipHFtj^il^ and those operators, where the top quark Pauli spinors are replaced by the 
antitop ones. As we show later in this section there are also contributions proportional to iTt to the 
photon and Z boson wave function renormalization constants, which have to be included to account 
for a proper treatment of the top quark decay phase space. We note that the velocity counting (|50|) 
implies that g Ug ^ v, where g is the SU(2) gauge coupling. Thus in a complete calculation of 
all electroweak effects one also has to determine 0{g'^) electroweak contributions and 0{a1) QCD 
corrections to the top decay width. In the framework of the Standard Model those corrections are 
known |^, 4C|. For our purposes, however, it is sufficient to consider the top quark width as an 
independent parameter, which we treat only to first order and which is not subject to higher-order 
corrections. 

The matching coefficients of the bilinear top spinor operators given above can be obtained 
by sandwiching the absorptive part of the top quark self energy Sj (in the full electroweak theory) 
between top quark Dirac spinors and expanding the result around the complex pole position. Keeping 
only the terms proportional to Tt, the result for the bilinear top Pauli spinor terms reads (k = (k^, k)) 



uik) 
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(51) 



The corresponding result for the bilinear antitop Pauli spinor terms reads 
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(52) 



where E^ = (fc^ + Mj?)!/^ and 
u{k) -- 
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(53) 



using the usual non-relativistic normalization for Dirac spinors. Thus in the presence of top decay we 
have to modify the NRQCD Lagrangian by adding the terms 



NRQCD 
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2Mf 



X- 



(54) 
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Physically, the dimension-5 operators multiplying the top width correspond to the time dilatation 
correction. This leads to the following modified versions of the Schrodinger equation ([l^ in momentum 
space representation 



AMf 
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iVt 1 



k' 



2 Mi 



G{k,k';q^ 



+ / -0^,V{k,p')G{p',k';q^) = (2vr)3 5(3)(fc-fcO 



(2 7r 



(55) 



and the free vertex functions {k) and G^ {k^ , k) in the integral equations (^5|), (^6|), (^0|) and (pl|): 
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(57) 



p2 

We note that it is not coercive to keep the term — ^ in the LO propagator, because it is of NNLO 
according to the power counting. We have adopted this convention because this choice leads to a 
simplification of the analytic form for the rest of the NNLO corrections in Eqs. (|56|) and (^). From 



the modified version of the Schrodinger equation (55) we can immediately see that the Fadin-Khoze 
replacement rule is valid at LO and NLO in the non-relativistic expansion, but it is inappropriate 
at NNLO. In fact it is not possible at all to consistently implement the top quark width at NNLO 
simply by shifting the centre-of-mass energy in a calculation for stable quarks. The resulting form of 
the optical theorem relations reads 
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(58) 



Im 



A'^iq^) 



2Nr 



4Nr 



Sk 

W) 



+ 00 



3 J 2tt 

— oo 



^|P(fe)|^r,. 



(2^)3 Ml 



(59) 



The first equality in Eqs. (58) and (|59| ) is an explicit integration over the phase space of the top decay 
products for off-shell top decay keeping only the term proportional to the top and antitop width for 



the top decay sub phase space. We emphasize again that in Eqs. (pSf ) and (|59|) no physical phase space 
boundaries for the integration over the top-antitop four momentum are implemented. All integrations 
are defined in the framework of the regularization scheme. 
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The reader might ask whether the inclusion of the top quark width into the top quark prop- 
agators leads to non-trivial modifications of the NRQCD short-distance coefficient C^.Q Because 
^ v"^ ^ a1, it is not inconceivable that there could be terms oc in C" at NNLO. To show that 
this is not the case we recall that C" is the modulus squared of the short-distance coefficient cf of 
the non-relativistic current t/j'^crx in Eq. (^). As mentioned in App. the short-distance coefficient 
Ci is the sum of amputated vertex diagrams in the full theory where the three momenta in the loop 
integrations are larger than the cutoff A ~ for ^/s = 2Mt. Thus, there is no integration over the 
top-antitop pole located at three momenta ~ po ~ Mtv, and we can conclude that the short-distance 
coefficient does not contain non-analytic terms involving the top quark width. Therefore can be 
expanded in iVt- For the same reason, the coefficients of an expansion in iVt do not contain a top- 
antitop cut and are real numbers. Because C" is the modulus squared of c\, the first non- vanishing 
term of an expansion of C" is proportional to {Tt/Mt)'^ ~ af, which is indeed beyond NNLO. How- 
ever, we emphasize that, in a complete calculation of all electroweak effects, C" will receive corrections 
Of 9^,g'^ ~ GFM^r, g' being the U(l) gauge coupling, which are formally of order Tt/Mt. These cor- 
rections do not come from the width contained in the top and antitop propagators but for instance 



from electroweak corrections to the vertices 7, Z ^ |41]. These corrections exist even in the case 
when the top quark is treated as a stable particle. 

If the top quark width is included in the NRQCD framework there is only one additional source 
of a linear dependence on Ff which comes from the phase space integration in Eqs. ( p8| ) and (59). In 



contrast to the case of a stable top quark, where the phase space is restricted to those top quark four 
momenta allowed by the centre-of-mass energy, the physical boundaries of the phase space integration 
/ (27r)^ ™ ""^^^ °^ unstable top quarks are determined by the allowed invariant masses of the top 
quark decay products. Assuming that the bottom quark is massless and that bottom quark and W 
boson are stable, and taking into account our routing convention (see Fig. |l|), the boundaries of the 
physical four-dimensional phase space integrations read (A;^ = 



d^k r d^k r dk'' 



(27r)4 J (27r)3 _ J 27r 

phase space 



(60) 



It is obvious that the physical limits of integration are not equivalent to the actual limits of integration 



on the RHS of Eqs. (58) and ( [59| ) as defined through our cutoff regularization scheme: 

(61) 
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For the total cross section on the LHS of Eqs. ( pq ) and (|59D, the difference between using the phase 
space integrations (|6^) or (|60|) can be expanded in and Fj, and can be accounted for by introducing 
additional photon and Z boson wave function renormalization constants into the NRQCD Lagrangian, 
which are proportional to iFf. The calculation of these counter-terms is straightforward. This leads 



^ One could equally well ask whether there is a need to introduce the operators j^ip'^crx ■^x^"'^ in the non- 
relativistic expansion of the currents, which produce and annihilate the top-antitop pair. At this point we prefer the 
language used in the text. 



19 



to the following form of the optical theorem relations for the vector- and axial-vector-current-induced 
correlators, 



Im 
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(63) 



where we have displayed the counter-terms for the phase space, each at the Born level, to first order 
in Tf and first order in the non-relativistic expansion. For simplicity we have set the bottom quark 
and the W boson masses to zero. As far as the velocity counting of the width is concerned, we in 
principle should have included also the 0(a<j) contributions for the counter-term in Eq. (|62[), because 
we formally have to count the vector-current-induced total cross section as order v. However, the 
counter-term contributions are only at the per cent level and not (yet) important phenomenologically, 
compared with the much larger uncertainties in the normalization of the total cross section from QCD, 



which are discussed in Sec. p.2| . We also note that the phase space counter-terms can be calculated 
entirely within the non-relativistic effective theory, and that there is no need to match to the full 
theory. This is because within the physical phase space boundaries the non-relativistic expansion of 
the phase space integration should be convergent.^ We also note that the calculation of the phase space 
counter-term for the S-P-wave interference, contributing e.g. in the top quark angular distribution, is 
in complete analogy to the calculation of the counter-terms in Eqs. ( |6^ ) and (|63|). We emphasize that 
the non-relativistic current correlators A", A"" and the vertex functions S and P in relations ( |6^ ) and 
(|6^) are still calculated in our cutoff scheme presented in Sec. ^. 

The expressions for the vector-current- and axial-vector-current-induced total cross sections 
valid at NNLO in the non-relativistic expansion, and properly including all effects of the top width at 
the Born level and leading order in the non-relativistic expansion, read 
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(64) 
(65) 



From Eqs. 
quarks {k 



|6^) and (^) we can derive the centre-of-mass three-momentum distributions of the top 
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(66) 



* This statement is a conjecture. In this context "convergence" for the phase space counter-terms is not meant to be 
associated with an expansion in Ft or po, but to the convergence in the coefBcient muhiplying the term hnear in Ft. AU 
higher-order terms in the non-relativistic expansion under the phase space integral can contribute to the term linear in 
Ft. 
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Unless \k\ is chosen close to the endpoint ((7^/4— M^)^/^ the numerical difference obtained by replacing 
the physical limits of the integrations by ±00 is negligible. We note that the three-momentum 
distributions shown in Eqs. ( |66D and ( |67[ ) are not equal to the physical observable three-momentum 
distributions, because the exchange of gluons between the top decay and production processes leads 
to additional non-negligible corrections at NLO |21, 22] and NNLO. These "interconnection" effects 
belong to the electroweak corrections, which are not treated in this work. The NNLO corrections 
to the three-momentum distribution calculated here are only a first step towards a complete NNLO 
treatment of the three- momentum distribution. We also want to mention that the three momentum 
distribution is strictly speaking an ambiguous quantity since the three (as well as the four) momentum 
of a coloured particle is an ambiguous concept. This is in contrast to the total cross section, which 
describes the rate of colour singlet top-antitop events. 

We conclude this section with some remarks on the inconsistencies that can arise if the Fadin- 
Khoze replacement rule "E' E+iTt" is employed at NNLO for the calculation of the top-antitop cross 
section close to threshold. We emphasize that there is nothing wrong, in principle, in calculating the 
current correlators for stable top quarks via the Schrodinger equation (II^), supplemented afterwards 



by the replacement 



2Mt 



— 2Mt + iVt- This corresponds essentially to the modification 
|n)(n| ^ 



En — E 



E-iVt 



(68) 



in the spectral representation of the Green function of Eq. (|T9| ) and is equivalent to keeping only the 
terms tpHTtxp and x^'^^tX in the modified version of the NRQCD Lagrangian. In this approach, also 
the optical theorem remains valid in the form 



Im 



(2vr)= 



l^(fc)rr, 



(69) 



for the vector current correlator, as an example. (For simplicity, we neglect the subtleties of the phase 
space effects, because they are irrelevant to this discussion.) However, there is a caveat, since it is 
possible, for the case of zero width, to simplify the form of Eq. ( |l9|) in a way that for the stable top quark 
case the results remain correct, whereas inconsistencies arise if the results undergo the replacement 
rule "-E — > E + iTf" ■ Such simplifications, based on the assumption that certain singular terms, which 
arise during the simplification, can be neglected, have in fact been carried out in Refs. [17, 11, 12, |l^. 
Li Refs. [17, 12, 34[ it was shown that the NNLO kinetic energy corrections and the Breit-Fermi 
potential in Eq. ([T9|), if they are treated as a perturbation to first order, can be rewritten in terms 
of an energy dependent Coulomb potential oc CpOis/Q^ x (pq/M^), a Darwin-like constant potential, 
and a potential tx ag/[Q[. Neglecting all NNLO corrections except the energy-dependent corrections 
to the Coulomb potential, the simplified version of the NNLO Schrodinger equation has the form|^ 



hi 

Mt 



Mt 



G{k,k') + 



v'''^{k-p')G{p\k') = {2^f 5^^\k - k') 



(2-kY 



(70) 



^ We emphasize that the neglect of the rest of the NNLO corrections does not affect the vahdity of the following 
arguments, because they are independent of the top quark width after the replacement rule '^E —> E + iVt" has been 
appUed. 
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where 



1 + 



3pg 



(71) 



Q2 V 2M2 

Equation ( [70[ ) is much easier to solve than the original Schrodinger equation ([T9|). For real energies, 
and if the corrections from the NNLO terms in Eqs. (^) and ( [70| ) are treated as a perturbation to 
first order only, the result obtained from Eqs. (^) and ( [70| ) are indeed equivalent, after a proper 
renormalization has been carried out. This was in fact the case for which the form of Eq. ( [70D has 
been derived in Refs. [0, 12, 34]. However, Eq. ( |70|) leads to inconsistencies for complex energies. 
This can be seen from the fact that for the total cross section calculated from Eq. (ff^), after applying 



the replacement rule j^^ ^ J4; + i^t , the actual form of the optical theorem relation reads 



Im^'' 



Mt 



(27r)3 J (27r 



\S{k)\^Tt - 6N, 



d^k r d^k' 



(27r)3 J (27r) 



S*{k)lm 



V'^'^ik - k') 



S{k'), (72) 



rather than Eq. (^). The additional term on the RHS of Eq. (|7^ ) originates from the energy dependent 
Coulomb- type potential in Eq. (|70|). However, the additional term does not correspond to any physical 
final state, because it corresponds to an absorption process in the potential. In other words, it is 
impossible to recover the total cross section from the momentum distribution, if one defines it as the 
integral over the physical final states represented by the first term on the RHS of Eq. (^). Including 
also the rest of the NNLO corrections not displayed in Eq. ([fO[), we have checked, with the numerical 
methods described in the next section, that the size of the second term on the RHS of Eq. ( [7^ ) is 

between about 5% (for - 2Mt ~ 5 GeV) and 20% (for - 2Mt 5 GeV) for the choices of 

parameters employed in the analysis of Sec. ^. (Similar results can already be obtained by analysing the 
known analytic solutions of the non-relativistic Coulomb problem |42, E^, E3] for a Coulomb potential 



with a complex coupling.) Thus for the determination of the total top-antitop cross section close to 
threshold, there is an unacceptable discrepancy between the integrated momentum distribution over 
physical final states and the absorptive part of the non-relativistic current correlator. We believe that 
the size of the second term on the RHS of Eq. (|7^ ) should in principle be taken as an estimate for the 
inherent uncertainties of using the simplified NNLO Schrodinger equation (|70|) supplemented by the 
replacement rule of Fadin and Khoze. We have checked, however, that the LHS of Eq. ([7^) is much 



closer to the correct result, obtained from the original Schrodinger equation (55), than the integrated 
momentum distribution. From this point of view the use of the simplified Schrodinger equation (70) 
might be justified for the total cross section, but is questionable for the momentum distribution. 



6 Numerical Implementation 



In this work we use numerical methods described in Refs. |^ 45, to determine the vertex functions 
S and P, which are the building blocks for the calculation of the total top-antitop production cross 
section and the three momentum distribution. Because the three- and four-dimensional versions of 
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the vertex functions are related through Eq. (42), it is sufficient to determine the amputated vertex 



functions S/G^ and P/G^ from the three-dimensional integral equations (35) and (^6|). The amputated 
vertex function are spherically symmetric and depend only on the modulus of the three momentum 



k. It is therefore possible to reduce Eqs. ( |35|) and (36) to one-dimensional integral equations 



Obviously, when solving Eqs. ( |35[ ) and (pq), the singular behaviour of the potentials V{k^p') 
and V^^°(fc,p') for p' ^ k requires special treatment. To avoid numerical problems we rewrite the 
integral equations for the amputated vertex functions as 



/C^(fc) 



S{k) 
Gf{k) 

1 + 7^72- I ^V{k,p')Gfip')IC%p') 



k" 


r <fp' 




1 (2^)3 




r (fp' 


6M2 . 


1 (27r)3 



IC'ik) 



^+6MF " ./ 7^mp')G^(p')(^^(p')-/C^(fc)) -B^ik))C''{k), (73) 
P(fc) 



G/(fc) 



1 - / ^V-^{k,p')Gfip') i^h^lC-{p')-IC'^{k)^-B^{k)!C^{k), (74) 



where B^ and B"" are defined as 

B'{k) ^ J^V{k,p')Gf{p'), (75) 

^■^(fc) ^ J ^V-°{k,p')Gf{p'). (76) 

As mentioned above, both /C^ and /C" depend only on the modulus of the three momentum. The 
angular dependence of the integrand in Eqs. (|73|) and (|7^ (including B^ and B"") is only coming 
from the potentials and the dot product kp' . The angular integration can be carried out analytically. 
The remaining one-dimensional integral equations are then solved numerically by discretization: the 
integrals / dp' {p' = \p'\) are transformed into sums over a fixed set of momenta p'*, and the 
integral equations for lC'"{k) and /C'^(fc) are each reduced to a system of linear equations, where the 
same set of momenta has to be used for the /c*.^ The resulting (complex) matrices are then inverted 
numerically to give the amputated vertex functions fC" and fC"" for the momenta A;'. 

In practice it turned out that the use of the Gaussian quadrature formulae is very efficient 
for the discretization. A surprisingly small number of points (of the order of 100) already leads to 
a high numerical accuracy. In addition, integrals were split into two (or more) parts and a suitable 
transformation of integration variables was applied wherever needed. It should also be noted that the 
finite width of the top quark is essential for the numerical stability of the method. It makes potentially 

^"At this point the subtraction carried out in Eqs. (^) and (^) becomes crucial. Without it the singularities in the 
potentials for p' k would be manifest for p^' = fc\ even in the case of integrable singularities. 
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dangerous denominators in the integrands of Eqs. ([73| ) and (^), which originate from the free Green 



function, weh behaved (compare Eqs. (37) and (l56|)). Clearly, the UV regularization by a momentum 
cutoff as discussed in Sec. Q is most naturally implemented in our numerical approach. It is sufficient 
to choose the momenta of the (Gauss-Legendre) grid to be limited by the value of the cutoff. Such a 
cutoff is, in principle, not needed in the case of pure Coulomb potentials. There the solution of the 
integral equations is possible without any cutoff and would correspond to a different regularization 
scheme with different short-distance coefficients. However, the potentials V^p, V^a and the kinematic 
corrections introduced at NNLO require a UV regularization already for purely numerical reasons, 
which can be seen from naive power counting in Eq. i^^. 



7 A First Analysis in the Pole Mass Scheme 

In this section we carry out a first brief analysis of the total cross sections Q^R" and R'^ and their 
three-momentum distributions in the pole mass scheme. We do this even though it is known that in 
the pole mass scheme there are uncomfortably large NNLO corrections in the location of the IS* peak 
position as well as in the normalization of the total vector-current-induced cross section K" [11, p^, [l^. 
We will show in Sec. ^ that these large corrections are a consequence of the pole mass scheme, and 
that the pole mass definition has to be abandoned as far as a precise extraction of a top quark mass 
from experimental data is concerned. Nevertheless, the pole mass is a well defined quantity in the 



framework of perturbation theory |4^, 48] , and, despite all its problems at larger orders of perturbation 
theory, remains a very convenient mass parameter to use for the formulation of (P)NRQCD and for 
calculations of the cross section. Thus a brief analysis in the pole mass scheme serves as a reference 
point with which results obtained with different approaches can be compared and from which we can 
visualize in which way alternative top quark mass definitions can improve the situation and in which 
way they cannot. 

In Figs. 1^ the total vector-current-induced cross section Q^R" is displayed for 344 GeV < 
a/^ < 352 GeV at LO (dotted hues), NLO (dashed lines) and NNLO (sohd lines). The LO curves 
are determined from Eq. (|55|), excluding all relativistic corrections and taking into account only the 
0{as) contribution to the Coulomb potential and setting = 1. In addition, the NLO curves also 
include the NLO corrections to and the 0{a'^) contributions of the Coulomb potential. At NNLO 
all contributions mentioned in this paper are included. Only at NNLO have we taken into account the 
phase space counter-terms displayed in Eqs. (|6^ ) and (l65|). In all figures shown in this section the top 
quark width is chosen as Tt = 1.43 GeV and the top quark pole mass as Mt = 175 GeV. Figure |3|a 
displays the dependence on the renormalization scale for /i = 15, 30 and 60 GeV for as{Mz) = 0.118, 
A = 175 GeV. For < 347 GeV, fj, = 15, 30 and 60 GeV corresponds to the upper, middle and 
lower curves, respectively. We note that in Refs. |11, 12, 13| the renormalization scale governing the 



strong coupling of the potentials in the Schrodinger equation (|19|) has been chosen between 50 and 
100 GeV. Regarding the fact that the inverse Bohr radius ~ MtUg is the scale that governs the cross 
section at NNLO, the natural renormalization scale is of order 30 GeV. This causes logarithms of the 
ratio A//i ~ l/a^ in the vector current correlator (see App. ^), which, however, are connected only 
to the running of the strong coupling. As has already been demonstrated in previous publications on 
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Figure 3: The total vector-current-induced cross section Q^R" for centre-of-mass energies 344 GeV < 
< 352 GeV in the pole mass scheme. The dependence on the renormalization scale /i (a), on the 
cutoff A (b) and on as{Mz) (c) is displayed. More details and the choice of parameters are given in 
the text. 



the NNLO corrections |11, 12, 13], the position of the 15 peak varies considerably at the different 
orders of approximation and for the different choices of the renormalization scale /x. (For explicit 
numbers see Table ||.) It is also evident that the normalization of the total cross section is subject 
to large corrections. We see that, in general, the NNLO corrections to the normalization are of order 
20% and as large as the NLO ones. Further, the dependence on the renormalization scale is even 
larger at NNLO than at NLO. In Fig. ^b the dependence of QfR" on the choice of the cutoff A 
is shown for as{Mz) = 0.118, ^ = 30 GeV and A = 90, 175 and 350 GeV. For LO and NNLO 
A = 90, 175 and 350 GeV correspond to the lower, middle and upper curves, respectively. For NLO 
A = 90, 175 and 350 GeV correspond to the lower, upper and middle curve, respectively. Whereas 
the dependence of the LO cross section on the choice of A is quite dramatic, because at LO there is 
no short-distance correction that could compensate for a variation in A, the variations at NLO and 
NNLO are significantly smaller (or order 5-10%). However, there is again no reduction of the variation 
from NLO to NNLO. The variation of the cross section with respect to A is small compared with the 
variation with respect to the renormalization scale fj, for centre-of-mass energies closer and below the 
peak. Figure lljc displays the dependence of Q^R" on the choice of as{Mz) for as{Mz) = 0.113,0.118 
and 0.123 and A = 175 GeV, /i = 30 GeV. As for the variations in the renormalization scale we see 
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Figure 4: The total axial-vector-current induced cross section i?" for centre-of-mass energies 
344 GeV < < 352 GeV in the pole mass scheme. The dependence on the renormalization scale 
H (a), on the cutoff A (b) and on as{Mz) (c) is displayed. More details and the choice of parameters 
are given in the text. 



that the position of the peak depends considerably on the choice of as- We observe a strong positive 
correlation between the choice of Mj and as{Mz)- Because the peak in the total cross section is the 
most pronounced feature of the total cross section, its behaviour directly reflects the quality of a top 
mass extraction from experimental data. Thus, if one would like to fit the pole mass to data for the 
cross section line-shape from a threshold scan, one finds a strong positive correlation between the pole 
mass and the strong coupling [Q] and a strong dependence of Mt on the choice of the renormalization 
scale leading to quite large theoretical uncertainties in the pole mass measurements. 

In Figs. |^,b,c the total axial-vector-current-induced cross section i?" is displayed for the same 
input parameters as in Figs. ^. Figure ^ shows the dependence on the renormalization scale for 
fj, = 15 (solid line), 30 (dashed line) and 60 GeV (dotted line). Figure ^ exhibits the dependence 
of the cutoff for A = 90 (solid line), 175 (dashed line) and 350 GeV (dotted line). Figure ^c shows 
ii" for as{M^) = 0.113 (solid line), 0.118 (dashed line) and 0.123 (dotted line). We observe that due 
to the V suppression of the axial-vector currents no peak-line enhancement as in the vector current 
case is visible. The variations of i?" with respect to the renormalization scale and the cutoff are quite 
large, because R*^ has only been determined at leading order, and the short-distance coefficient 
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Figure 5: The three- momentum distribution of the vector-current-induced cross section QfR" for 
centre-of-mass energies = Mp^ak and Mp^ak + 5 GeV in the pole mass scheme. The dependence 
on the renormahzation scale /i (a), on the cutoff A (b) and on as{Mz) (c) is displayed. More details 
and the choice of parameters are given in the text. 



does not contain any corrections that could compensate for the variations. We note that the NLO 
corrections to i?" can be implemented in the same way as for K" . They are, however, beyond NNLO 
and have therefore not been included into our analysis. From the phenomenological point of view the 
next-to-leading order corrections are irrelevant if one takes into account the small normalization of -R" 
compared to i?". 

In Figs. |5|a,b and c the LO (dotted lines), NLO (dashed line) and NNLO (solid lines) top- 
antitop vector-current-induced three momentum distribution Q^dR" /d\kt\ is shown for < \kt\ < 
60 GeV and both for centre-of-mass energies exactly on top of the visible peak, ^/(p = Mp^ak , and for 
y/q^ = Mpeak + 5 GeV. The input parameters have been chosen as in Figs. ^. Figure ||a shows the 
distributions for = 15 and 60 GeV. At LO and NNLO = 15 GeV corresponds to the upper curves 
below the peak and /x = 60 GeV to the lower curves. At NLO /i = 60 GeV corresponds to the higher 
peak and fi = 15 GeV to the lower. Figure |5|b displays the dependence of the distributions on the 
cutoff for A = 90 (lower curves) and 350 GeV (upper curves), and Fig. |5|c exhibits the dependence of 
the distributions on the strong coupling for as{Mz) = 0.113 and 0.123. Below the peak, the larger 
value of as{Mz) always corresponds to the upper curve. As for the total cross section, we observe 
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Figure 6: The three-momentum distribution of the axial-vector-current-induced cross section 72" for 
centre-of-mass energies = Mpeak and Mp^ak + 5 GeV in the pole mass scheme. The dependence 
on the renormalization scale /i (a), on the cutoff A (b) and on as{Mz) (c) is displayed. More details 
and the choice of parameters are given in the text. 



a strong dependence of the normalization of the distributions on the renormalization scale and the 
cutoff. For = Mpeafc + 5 GeV, the dependence of the peak position on the renormalization scale 
is particularly strong]^ 

Finally, in Figs. |6|a,b and c the top-antitop axial- vector-current-induced three momentum dis- 
tribution dR"" /d\kt\ is shown for < \kt\ < 60 GeV and for centre-of-mass energies exactly on top of 
the visible peak, ^/(p = Mp^ak, and for = Mpeak + ^ GeV. The input parameters have been chosen 
as before. Figure ^ shows the distributions for fi = 15 (solid curves), 30 (dashed curves) and 60 GeV 
(dotted curves) . Figure ^b displays the dependence of the distributions on the cutoff for A = 90 (solid 
curves) 175 (dashed curves) and 350 GeV (dotted curves), and Fig. |^ exhibits the dependence of the 
distributions on the strong coupling for as{Mz) = 0.113 (solid curves), 0.118 (dashed curves) and 
0.123 (dotted curves). As expected, the momentum distribution is strongly suppressed for smaller 
centre-of-mass energies. The variations of the normalization of the distributions are comparable to 
the variations of the normalization of the total cross sections. 

The peak position is always located approximately at |fet| ~ (pg + )^'^'*, which can be regarded as the effective 

three momentum of the top quarks. 
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8 Theoretical Uncertainties 



In the previous section we have seen that the location of the peak position and the normaHzation of 
the vector-current-induced total cross section, as well as of its three-momentum distribution, receive 
large NNLO corrections in the pole mass scheme. In the following two subsections we discuss the 
origin of the large corrections in the vector-current-induced total cross section. For the peak position 
of the total cross section we propose a solution, which allows for a considerable stabilization. 



8.1 The IS Top Quark Mass and the Peak Position 

Prom past experience in the theoretical description of i?-meson decays [EM 5C], it is well known that 



the pole mass, defined as the pole of the perturbative quark propagator, although infrared-finite to 



all orders in perturbation theory ||47|, |4^, is a concept that is ambiguous to an amount of order 
^QCD- This might be a reflection of the fact that the perturbative quark pole does not exist in reality 
because of confinement. Within the framework of perturbation theory this ambiguity is caused by an 
n-factorial increase of the coefficients in the perturbative relation between the pole mass and a short- 
distance mass such as MS. The large corrections are caused by an increasing infrared sensitivity of 
the perturbative coefficients for large orders. It has also been shown that the large top width does not 
lead to a suppression of these large corrections From this point of view the unstable behaviour 
of the peak position in the total cross section is not unexpected and it would be quite appealing 
conceptually to conclude that the use of a short-distance mass instead of the pole definition would 
cure the problem. 

In fact, it has been demonstrated in Refs. ||5^, |5^ that the dominant source of infrared sen- 
sitivity in the Green function of the Schrodinger equation ( |l9|) comes from the terms in the static 



energy 2Mt + V{r). Whereas the rest (pole) mass energy 2Mt and the potential energy V (which has 



traditionally been calculated in the pole mass scheme |26, 27, |2^) are individually ambiguous to 
an amount of order Aqcd |5^> H], the sum of both is not ||52|, |5^. This shows that quantities 
such as spectra or the total cross section calculated from Eq. (|l^) are much less sensitive to infrared 
momenta than the pole mass itself, rendering it an irrelevant mass parameter. Thus, any sensible 
mass definition used to parametrize the top-antitop cross section close to threshold should have no 
ambiguity of order Aqcd, i-e. it should be a short-distance mass. 

However, we emphasize that, in practice, large corrections at lower orders in perturbation 
theory in the pole scheme do not necessarily come from the ambiguity in the pole mass. This is because 
the cancellations of infrared sensitive contributions in static and rest mass energy is a phenomenon 
that is relevant in large orders where the corresponding series are dominated by the most infrared 
sensitive contributions in the loop integrals. Thus, large corrections in perturbation theory at low 
orders could very well come from scales which are not infrared. To get a clearer picture for the case of 
the peak position in the total cross section, let us have a look at the size of the individual corrections 
to the peak position. In Table |l] we have displayed the LO, NLO and NNLO contributions to the peak 
position with respect to the pole rest mass energy 
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Table 1: LO, NLO and NNLO contributions to the peak position of the total vector current induced 
cross section R" in GeV in the pole mass scheme for Mt = 175 GeV, Tt = 1.43 GeV, as{Mz) = 
0.113,0.118 and 0.123, and fi = 15,30 and 60 GeV, respectively. For the strong coupling two-loop 
running has been employed. The results are insensitive to the choice of the cutoff scale A ~ Mf. 



= 2Mt - 5Mp,ak , (77) 

for Mt = 175 GeV, Tt = 1.43 GeV, as{Mz) = 0.113,0.118 and 0.123, and = 15,30 and 60 GeV. At 
NNLO we have separated from the rest the contributions with the highest power of /3o, which represent 
the contributions most sensitive to infrared momenta^. All the large-/?o terms originate from the 
Coulomb potential 14, Eq. (p2|). From the numbers presented in Table |l| we see that depending on 
the choice of the renormalization scale the large-/3o contributions at NNLO contribute between about 
30 and 60% to the total NNLO corrections to the peak position. Thus at NNLO the large shift in the 
peak position consists to approximately equal parts of corrections very sensitive to infrared momenta 
and corrections coming from subleading infrared terms and relativistic corrections. From this we can 
conclude that using some unspecified short-distance mass definition instead of the pole mass does 
not necessarily lead to smaller NNLO corrections to the peak position because the most infrared 
sensitive terms are not yet dominating at NNLO. For the same reason, we cannot conclude that some 
unspecified short-distance mass definition necessarily leads to a significantly reduced renormalization 
scale dependence of the NNLO peak position or to a smaller correlation between the peak position 
and the strong coupling. Thus, the question of which mass definition one should use is not only a 
conceptual issue, but also a practical one. 

We formulate two requirements for a proper top mass definition for the total cross section close 
to threshold: 

A) it must be a short-distance mass, and 

These contributions are determined by using the replacement rule n/ — + — |/3o for the terms with the highest power 
of 71/, where is the number of light quark species. This method is called "naive non-Abelianization" and accounts for 
the most infrared sensitive contributions in perturbative series relating the pole mass to a short-distance mass . The 
result of this replacement is referred to as the "large-/3o" limit. 
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B) it must lead to a considerable stabilization of the peak location with respect to the order of 
approximation used and also to variations of parameters such as the strong coupling or the 
renormalization scale. 



Requirement A reflects the necessity that, if a top mass determination at the linear collider with 
uncertainties of 200 MeV or even better is intended, the corresponding mass parameter must be free 
of intrinsic ambiguities of order Kqcd- In addition, only a short-distance mass can be reliably related 
to the MS top quark mass, which is the preferred mass parameter used in calculations at high energies 
and for top quark corrections to electroweak precision observables. Requirement B ensures that the 
mass parameter can be extracted from experimental data with small systematic uncertainties. 

The mass that seems to be most appropriate to us to fulfil this task, because it is closely 
related to the peak position in the vector-current-induced cross section, is what we call the IS mass. 
Mis- The IS mass is defined as half the perturbative mass of the fictitious toponium l^Si ground 
state, where the top quark is assumed to be stable. Expressed in terms of the pole mass, the IS mass 
reads (a^ = as(/^)) fH, 57 1: 
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where 



L = In 



Cf as Mt 



(78) 



(79) 



and the contributions at LO, NLO and NNLO are labelled by the powers e, and e'^, respectively, 
of the auxiliary parameter e = 1. In general, the electroweak corrections not calculated in this work 
can lead to further corrections in Eq. (|78|). We note that 2Mis is not equal to the actual peak 
position visible in the vector-current-induced total cross section because the top quark width leads to 
an additional shift of the peak by about +200 MeVj^ In principle it would also be possible to define 
a mass that would be equal to half the actual visible peak position. Except for additional corrections 
coming from the top quark width, this would also require the inclusion of an additional shift coming 
from the axial-vector-induced cross section. Such a definition would, however, not necessarily be more 
useful, since the experimentally measurable line-shape of the total cross section at a future e~^e~ or 
fj-^fj-^ collider will be distorted by initial state radiation and beamstrahlung. In the case of a muon 
collider these effects lead to an additional shift and in the case of the e~^e~ linear collider even to 
the disappearance of the peak P]. Thus one has to consider the 15 mass, like the MS mass, as a 



fictitious mass parameter, which to NNLO is defined through the perturbative series given in Eq. (78). 



The difference between 2Mis and the peak location of the vector-current-induced total cross section is proportional 
to Ft times a function of Tt/iMta'i). For Ft < Mta^ the difference Mpeak - 2Mis is proportional to /{MtCpal)'^ . 
The size of the difference between 2Mis and the peak location of about 4-200 MeV can be seen in Figs. ^. 
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Nevertheless, twice the IS mass is quite close to the peak location and, by construction, leads to a 
considerable reduction of the variation of the peak position with respect to the order of approximation, 
the strong coupling and the renormalization scale. 

To show that Mis is indeed a short-distance mass we recall that the static energy 2Mt + 
14 (r) represents the dominant infrared sensitive contribution in the Schrodinger equation (|T9[). The 
difference between Mis and the pole mass isf^ 

Mis-M, = I J -^^^^lsiP)n{p,q)^,s{q), (80) 

where TC is the Hamiltonian of Eq. ( [T9| ) and <I>is the normalized wave function of the IS* state in 
momentum space representation. The dominant infrared sensitive contribution in relation (|80|) reads 



{Mis-M,)''' ~ \J-§^^^*s{Pmp-q)^ls{q). (81) 



Because the infrared region in Eq. (|8lD is given by |p — q| < /xj, where /x/ is much smaller than the 
inverse Bohr radius ~ MtUg, the characteristic scale governing the dynamics described by the wave 
function, we can simplify relation (^T]), 



[Mis 



It has been shown in Refs. |5^, |5^ that the RHS of Eq. (^) is equivalent to the dominant infrared 
contributions of the difference between the MS and pole mass. Therefore the relation between Mis 
and the MS mass ifit only contains subleading infrared contributions, which are suppressed by at least 
one power of l/M^. In other words the ambiguity in the relation between Mis a-nd rht is parametrically 
of order K^(j^/Mt. This proves that Mis is a short-distance mass. We also see from Eq. ( p2[ ) that, if 
the pole mass is expressed in terms of the 15" mass and if the resulting mass difference 2(Mf — Mis) is 
absorbed into the potential, the rest mass and the potential energy term contained in the total static 
energy are individually free of ambiguities of order Aqcd- The RHS of Eq. (|8^ ) just subtracts the low 
momentum (i.e. dominant infrared sensitive) contribution from the Coulomb potential Vc{x). 

We note that in order to implement the IS mass definition into our numerical codes, which 
solve the integral equations ( |35|) and (^6|), we have to invert relation (^8|). It has been shown in 
Refs. pO| , |6l| that a consistent way to achieve this task is to carry out the inversion with respect to 
the auxiliary parameter e. For the reason that this modified perturbative expansion has been applied 

Strictly speaking, the simple form of Eq. ( po| ) is true only up to NNLO because of retardation effects, which set in 
at N'^LO. Thus, in general, there would also be a non-trivial integration over time components. The form of our proof 
also depends on the assumption that the static potential is an infrared finite quantity. That this is most probably not 
the case was already pointed out some time ago in Ref. [|8|, because the perturbative static potential might become 
sensitive to scales below the inverse Bohr radius at ©(qj). Some contributions at 0(Qs) have recently been calculated 
in Ref. [^. Up to 0{a^), i.e. NNLO in the non-relativistic expansion, the perturbative potential has been proven to 
be finite by complete calculations Because Mis is defined as a physical quantity this would not affect the final 

conclusion that it is a short-distance mass, but it would change the form of the proof considerably. 
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Figure 7: The total vector-current-induced cross section Q^R" for centre-of-mass energies 346 GeV < 
< 354 GeV in the 15 mass scheme. The dependence on the renormalization scale ^ (a), on the 
cutoff A (b) and on as{Mz) (c) is displayed. More details and the choice of parameters are given in 
the text. 



for the first time to express inclusive B decays in terms of the T(15') mass, it has been called the 
"Upsilon expansion". ^] If the 15 mass is expressed in terms of the MS mass, which is related 
to the pole mass by a series of the form rht — Mt = Mt J2'^=i ^n"" > one has to consider a term oc a" 
in this relation of order e" in the Upsilon expansion. In other words, if one relates the 15 mass to a 
mass which is different from the pole mass, one must combine terms of different order in Og. As an 
example, this means that in order to relate the N'^LO 15 mass to the MS mass one needs to know its 
relation to the pole mass to 0{ag~^^). This is necessary because this is the only way in which the high 
order large perturbative corrections coming from infrared-sensitive terms are cancelled. 

In Figs. 1^ the total vector-current-induced cross section Q^R" is displayed in the 15 scheme 
for 346 GeV < < 354 GeV at LO (dotted lines), NLO (dashed hues) and NNLO (sohd lines). In 
all figures shown in this section the top quark width is chosen as F^ = 1.43 GeV and the top quark 15 
mass as Mis = 175 GeV. Figure ^ displays the dependence on the renormalization scale for fi = 15, 
30 and 60 GeV for a,(Mz) = 0.118 and A = 175 GeV. At LO and NNLO the choices n = 15, 30 
and 60 GeV correspond to the upper, middle and lower curves. At NLO the choices /i = 15, 30 and 
60 GeV correspond to the lower, middle and upper curves for centre-of-mass energies below the peak 
position. In Fig. |^ the dependence of Q^R" on the choice of the cutoff A is shown for as{Mz) = 0.118, 
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Figure 8: The total axial-vector-current-induced cross section i?" for centre-of-mass energies 
346 GeV < \fcp' < 354 GeV in the IS* mass scheme. The dependence on the renormalization scale fi 
(a), on the cutoff A (b) and on as{Mz) (c) is displayed. More details and the choice of parameters 
are given in the text. 

= 30 GeV and A = 90 (lower curves), 175 (middle curves) and 350 GeV (upper curves). Figure ^c 
displays the dependence of Q^R" on the choice of as{Mz) for as{Mz) = 0.113 (lower curves), 0.118 
(middle curves) and 0.123 (upper curves) and A = 175 GeV, /i = 30 GeV. Comparing the result with 
the curves displayed in Figs. ^, the improvement of the stability of the peak position is evident. The 
strong dependence on the renormalization scale and the strong correlation with as{Mz) have vanished. 
However, we also observe that the large corrections in the normalization of the curves are essentially 
not affected at all. 

In Figs. ^,b,c the total axial-vector-current-induced cross section R"' is displayed in the IS 
mass scheme for the same input parameters as in Figs. ^. Figure ^ shows the dependence on the 
renormalization scale for = 15 (solid line), 30 (dashed line) and 60 GeV (dotted line), respectively. 
Figure |8|b exhibits the dependence of the cutoff for A = 90 (solid line), 175 (dashed line) and 350 GeV 
(dotted line). Figure |c shows for as(M^) = 0.113 (solid line), 0.118 (dashed line) and 0.123 
(dotted line). Compared to the plots in the pole mass scheme, we find a slightly smaller variation in 
the normalization with respect to the renormalization scale and the choice of as{Mz)- Clearly the 
effects of using the IS scheme instead of the pole one are much smaller in the axial- vector case because 
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Figure 9: The three-momentum distribution of the vector-current-induced cross section QfR" for 
centre-of-mass energies = Mp^ak and Mp^ak + 5 GeV in the IS* mass scheme. The dependence on 
the renormahzation scale (j, (a), on the cutoff A (b) and on as{Mz) (c) is displayed. More details and 
the choice of parameters are given in the text. 



no peak is visible there. 

In Figs. Pa,b and c the LO (dotted lines), NLO (dashed lines) and NNLO (solid lines) top- 
antitop vector-current-induced three-momentum distribution Q^dK" / d\kt\ is shown for < \kt\ < 
60 GeV in the 15 mass scheme for centre-of-mass energies exactly on top of the visible peak, ^/cp = 
Mpeak and for ^/q^ = Mp^ak + 5 GeV. The input parameters have been chosen as in Figs. 0. Figure |9|a 
shows the distributions for ^ = 15 and 60 GeV. At LO and NNLO /i = 15 GeV corresponds to the 
upper curves and n = 60 GeV to the lower curves for centre-of-mass energies below the peak. At NLO 
// = 60 GeV corresponds to the higher peak and = 15 GeV to the lower. Figure |9|b displays the 
dependence of the distributions on the cutoff for A = 90 (lower curves) and 350 GeV (upper curves), 
and Fig. ^ exhibits the dependence of the distributions on the strong coupling for as{Mz) = 0.113 
and 0.123. Below the peak the larger value of as{Mz) always corresponds to the upper curve. 

In Figs. 10i,b and c the top-antitop axial- vector-current-induced three-momentum distribution 
dR°'/d\kt\ is shown in the IS scheme for < |fcf| < 60 GeV and both for centre-of-mass energies exactly 
on top of the visible peak, ^/(p = Mpeak, and for ^/(p = Mpeak + 5 GeV. The input parameters have 
been chosen as before. Figure ICa shows the distribution for = 15 (solid curves), 30 (dashed curves) 
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Figure 10: The three-momentum distribution of the axial-vector-current-induced cross section i?" for 
centre-of-mass energies = Mp^ak and Mp^ak + 5 GeV in the 15* mass scheme. The dependence on 
the renormahzation scale (j, (a), on the cutoff A (b) and on as{Mz) (c) is displayed. More details and 
the choice of parameters are given in the text. 



and 60 GeV (dotted curves). Figure |TO|b displays the dependence of the distribution on the cutoff for 
A = 90 (solid curves), 175 (dashed curves) and 350 GeV (dotted curves), and Fig. px| c exhibits the 
dependence of the distribution on the strong coupling for as{Mz) = 0.113 (solid curves) 0.118 (dashed 
curves) and 0.123 (dotted curves). The curves shown in Figs. 10 are somewhat higher than in Figs. 
because the choice of 175 GeV for the top quark mass corresponds to a higher value for Mp^ak in the 
IS scheme. From Figs. |9| and [l^ it is evident that the IS" scheme does not essentially affect at all the 
three-momentum distributions. Compared to the results in the pole mass scheme the variations of 
the peak position remain unchanged. This can be understood from the fact that a mass redefinition 
corresponds to a shift in the centre-of-mass energy, but leaves the definition of the off-shell top quark 
three-momentum unchanged. 

In Table |2| we have displayed the LO, NLO and NNLO corrections to the peak position with 
respect to 2M15: 

Mpeak = 2M,s + 5M^,%^s + ^M^eak,lS + ^M^eakfs 

= 2Mis + 6Mp,ak,is , (83) 
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/i[GeV] 


as{Mz) 










15 


0.113 


0.21 


0.03 


-0.03 


0.20 


30 




0.38 


-0.11 


-0.09 


0.17 


60 




0.78 


-0.50 


-0.11 


0.17 


15 


0.118 


0.16 


0.02 


-0.00 


0.17 


30 




0.30 


-0.09 


-0.08 


0.12 


60 




0.54 


-0.32 


-0.10 


0.11 


15 


0.123 


0.12 


-0.00 


0.04 


0.16 


30 




0.23 


-0.07 


-0.08 


0.08 


60 




0.42 


-0.26 


-0.10 


0.07 



Table 2: LO, NLO and NNLO contributions to the peak position of the total vector-current-induced 
cross section K" in GeV in the 15 mass scheme for Mis = 175 GeV, Tt = 1.43 GeV, as{Mz) = 
0.113,0.118 and 0.123, and fi = 15,30 and 60 GeV, respectively. For the strong coupling two-loop 
running has been employed. The results are insensitive to the choice of the cutoff scale A ~ 175 GeV. 



in the 15 mass scheme for Mis = 175 GeV, Ft = 1.43 GeV, as{Mz) = 0.113,0.118 and 0.123, and 
^j, = 15, 30 and 60 GeV for various choices of the renormalization scale fi and the strong coupling 
as{Mz)- Taking the size of the NLO and NNLO corrections as a measure for the present theoretical 
uncertainty in the peak position, and assuming that the latter can be used to estimate the theoretical 
uncertainty in the determination of Mis, we find that this uncertainty is approximately 200 MeV. 
If the effects of beamstrahlung and initial state radiation at a future e~*'e~ or muon pair collider do 
not spoil a precise determination of the 15 mass with an uncertainty of 200 MeV one has to ask 
the question how Mis is related to the top mass parameters usually used for calculations of physical 
observables that are not related to the threshold regime. In principle, one could use Mis as a new top 
mass parameter in its own right. This would, of course, require that all formulae be expressed in terms 
of Mis, using the Upsilon expansion discussed after Eq. (^2|). A more economical way is to relate the 
15 mass to the MS top quark mass, which is a common mass parameter for perturbative calculations 
involving heavy quarks and which, in a number of cases, even leads to improved convergence properties 
of the perturbative series.|^ Because the 15 mass is a short-distance mass, its perturbative relation 
to the MS mass is much better behaved at large orders than the corresponding relation of the pole 
mass. The relation between Mis and fht{fht) can be derived from Eq. (|7^) and the relation between 
pole mass and fht{rht) using the Upsilon expansion discussed above. We emphasize that the three- 
loop relation between the pole and the MS mass is needed to relate M15 and rht{rht) at NNLO 
accuracy. Assuming that those 3-loop corrections can be approximated by the known corrections in 
the large-/3o limit we find the following numerical value for ?fif (mj ) for Mis = 175 it 0.2 GeV and 

^® Prominent cases are the top quark QCD corrections to the p-parameter and the massive quark pair production 
cross section at large energies 
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as{Mz) = 0.118 ibxO.OOl, 



mt{mt) 



175 - 7.58 e(LO) - 0.96e2(NLO) - 0.23e^(NNLO) 



±0.2(5Mi5) ±x0.07(fc. 



GeV. 



(84) 



For the numbers given in Eq. ( p4D we have assumed an uncertainty in the value of the strong couphng 
at Mz of X 0.001 in order to demonstrate the importance of as for the determination of 'mt{rht). This 
uncertainty is independent of the order to which the relation between the pole and the MS mass is 
known because it comes from the LO term. We note that this fact shows that the strong correlation of 
the peak position to the strong coupling, which was visible in the pole mass scheme, is not necessarily 
eliminated by adopting the 15 scheme. This correlation might come back whenever the IS" mass is 
related to another short-distance mass or is used as a parameter in other quantities. However, the use 
of the 15 mass has the advantage to free the process of the mass extraction from the total cross section 
close to threshold also from strong dependences on other parameters such as the renormalization scale 
or the order of approximation used. Therefore systematic uncertainties are expected to be smaller if 



the 15 scheme is used for the threshold calculations. Equation ( p4D shows that the knowledge of the 
3-loop corrections in the relation of pole and MS mass and a small uncertainty in as{Mz) are crucial 
for a determination of fht{fht) with uncertainties comparable to SMis- 

In recent literature there have been two other proposals for alternative short-distance mass 
definitions, which can also be used for a measurement of the top quark mass from the total cross 
section. In Refs. |64, ^ the "low scale running mass" was proposed to subtract the infrared behaviour 
from the heavy quark self energy. The "low scale running mass" was devised in order to improve the 
convergence of the perturbative series describing the contributions leading in l/Mb in inclusive B- 
meson decays. Due to the universality of the dominant infrared sensitive contribution, the low scale 
running mass can also serve as a top mass definition, which leads to an improved stability of the peak 
position in the total cross section. The low scale running mass depends on the cutoff ^ls^ which limits 
the momenta that are subtracted from the self energy. At order (i.e. at LO) its relation to the 



pole mass reads [64 



mP(/iL5) 



Mt 



16 as 

-T. ^L5 

9 vr 



l + 0(a,) + 0( 



(85) 



By adjusting the scale ^ls in such a way that the RHS of Eq. ( |85| ) is comparable in size to the 
RHS of Eq. ( |78D the position of the peak in the total cross section can be stabilized. In Ref. |53] the 
"potential-subtracted" mass was proposed. It subtracts the dominant infrared-sensitive contribution 
in the Schrodinger equation (|l9|), which is contained in the static potential Vc- The subtraction is in 
fact equal to the RHS of Eq. (|8^). Like the low scale running mass, the potential-subtracted mass 
depends on a cutoff, ups- At order as (LO) the relation to the pole mass reads 



mP(^ 



PS) 



Mt 



4 as 
3 vr 



l + 0(as 



(86) 



As for the low scale running mass, the scale fips can be adjusted in such a way that the RHS of Eq. (86) 
is comparable in size to the RHS of Eq. (p^). To achieve this, fj,ps has to be chosen of the order of the 
inverse Bohr radius ~ Mjag, which is much larger than the scale /i/ *C Mtas introduced in Eq. (|82|) . 
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For ^ps = 2i^LS the low scale running and the potential-subtracted mass lead to approximately 
equivalent results. However, the stabilization of the peak position can be expected to be slightly worse 
than for the 15 mass if ^ps or l/iLS are not fine-tuned. In addition, the results that could finally 
be obtained for the MS top mass can depend on the value that is chosen for the cutoff scale ^ls and 
^J'PS■ 



8.2 Normalization of the Total Cross Section 

In the previous subsection we have demonstrated that a proper redefinition of the top quark mass 
leads to a considerable improvement in the stability of the peak position in the vector-current-induced 
total cross section R" . However, there have been only marginal changes in the size of the NNLO 
corrections to the overall normalization of the line-shape. Compared to the NLO normalization of 
the total vector-current-induced cross section, the NNLO corrections are between 15 and 25%, which 
is rather large if one recalls that the NNLO corrections are parametrically of order ~ a^. In 
this section we try to find some answers to the question, whether the large NNLO corrections to 
the normalization of the total vector-current-induced cross section have to be interpreted as a sign 
that the non-relativistic expansion for the top-antitop cross section close to threshold breaks down. 
Clearly, this question can only be answered reliably after the complete N^LO corrections have been 
determined, which are, unfortunately, beyond the capabilities of present technology. We therefore 
analyse the NNLO corrections to the normalization of the total cross section with respect to their 
sensitivity to infrared momenta and carry out a comparison to the one- and two-loop cross section for 
energies far above the top-antitop threshold, where conventional perturbation theory in is believed 
to be reliable. We provide arguments that the large NNLO corrections to the normalization are genuine 
0{v^,a1) relativistic corrections, which cannot be removed by changing the definition of or the 
top quark mass, and that their size does not necessarily indicate a breakdown of the non-relativistic 
expansion used in this work. 

As far as a redefinition of the top quark mass is concerned, it is quite obvious that it cannot sig- 
nificantly affect the normalization of the total cross section because the dominant effect in a mass shift 
is an energy shift of the entire line-shape. Nevertheless, it is quite interesting that the normalization 
is at all insensitive to the dominant infrared-sensitive terms in the Schrodinger equation (|T^)P^, which, 
in the pole mass scheme, would cause the corrections to the peak position to grow factorially at large 
orders of perturbation theory. To show this let us recall that the total vector-current-induced cross 
section is proportional to the absorptive part of the Green function, with both arguments evaluated 
at the origin in configuration space representation: 



where the sum extends over discrete and continuum states with S wave quantum numbers. Thus, 
for fixed energy the normalization only depends on the wave function. Repeating the steps following 



In our case the two issues are in fact connected to each other. But it is important to conceptually separate the issue 
of a simple energy shift from the more fundamental question of infrared sensitivity. 
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Figure 11: The total vector-current-induced cross section QfR" close to threshold for Mis = A = 
175 GeV, /i = 30 GeV, OsiMz) = 0.118 and Tt = 1.43 GeV at NLO (dotted curve) and NNLO (solid 
curve). The dash-dotted curve is NLO including also the NNLO corrections to the Coulomb potential 
Vc, and the dashed line contains, in addition, all Abelian NNLO corrections. The differences between 
the curves indicates the size of individual NNLO relativistic corrections. 



Eq. (80) we find that the correction to the wave function coming from the dominant infrared-sensitive 
terms in the Schrodinger equation reads 



<5$„(0) 
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0, 



i.e. it vanishes because of the orthogonality of the wave functions. Therefore the large corrections in 
the normalization of the total cross section are not related to an infrared sensitivity of the corrections, 
in particular at large orders. To demonstrate that this is also the case for the NNLO corrections 
calculated in this work, we have displayed in Fig. ^ the total vector-current-induced cross section 
QlR" close to threshold for Mis = A = 175 GeV, /i = 30 GeV, as{Mz) = 0.118 and Tt = 1.43 GeV 
successively including various NNLO corrections. The dotted line represents the NLO cross section and 
the solid line the NNLO one. The dash-dotted line is the NLO cross section including also the NNLO 
corrections to the Coulomb potential Vc, the dashed line contains, in addition, all Abelian NNLO 
corrections, i.e. those that do not involve the SU(3) group theoretical factor Ca- The separation 
of the NNLO corrections into those coming from the Coulomb potential and from Abelian and non- 
Abelian relativistic corrections is gauge- invariant. The difference between the dashed and the solid 
curve represents the corrections of the non-Abelian NNLO effects originating from the potential V^a 
and those 0{a1) contributions to the short-distance coefficient that are proportional to Ca- From 
the rather small difference between the dotted and the dash-dotted curves (2-4%) we see that the large 
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NNLO corrections to the normalization are not related to the corrections in the Coulomb potential. 
Because a redefinition of the strong coupling would mainly affect the size of the higher-order corrections 
in the Coulomb potential, we can conclude that using a different scheme for the strong coupling (such 
as the ^-scheme []66| , |67| ) will not significantly affect the size of the NNLO corrections. The curves 
plotted in Fig. 11 demonstrate that the 0(20%) NNLO correction to the normalization is a sum of 
corrections, each of which positive and individually either smaller than or approximately equal to 
0(10%). Although this observation, of course, cannot be taken as a proof that the still unknown 
N^LO corrections are smaller than the NNLO ones, it indicates that the size of the latter does not 
necessarily have to be taken as an argument for the non-relativistic expansion to break down for the 
normalization of the total cross section. 

An interesting insight into the question of how to interpret the large normalization corrections 
can also be obtained by comparing the total cross section line-shape, which we have calculated in 
the threshold regime, with earlier calculations of the total cross section for higher energies, where a 
resummation of Coulomb singular terms is not yet necessary and perturbation theory in is believed 
to be reliable. ||6^ We would like to note that it is the large mass of the top quark that allows us 
to draw conclusions from a comparison of the threshold cross section with the one calculated for 
higher energies. To illustrate this we recall that our calculation of the threshold cross section is valid 
if the hierarchy as,v <C 1 is satisfied, where the scale of the strong coupling is of the order of the 
inverse Bohr radius, the kinetic energy, or the top width. This means that the threshold cross section 
represents a simultaneous expansion in a<j and where powers of (a<j/f) are resummed to all orders 
in Us- The high energy cross section, on the other hand, is valid if <C v, 1, where the scale in the 
strong coupling if of order the top-antitop relative momentum or the centre-of-mass energy. Thus a 
comparison of the threshold results with the high energy perturbative ones is only sensible if there 
exists a kinematic regime where both hierarchies are satisfied at the same time, i.e. if <C w <C 1. 
In this regime the effects of the resummation of powers of (a^/f) not contained in the high energy 
cross section should be small as well as the effects of velocity corrections beyond NNLO, which are not 
contained in the threshold cross section. Obviously this relation is difficult or impossible to satisfy for 
bottom or charm quarks, but it is possible for the top quark case. For as{Mtas) ~ 0.13 we can argue 
that a meaningful comparison between threshold and high energy cross section should be possible 
for V ~ 0.3-0.4, which corresponds to ^/s « 365 GeV. In Fig. |l2| we have plotted the threshold and 
the high energy cross sections at NLO/NNLO and 0{as) /0{a1), respectively, for the renormalization 
scales /i = 25 GeV (dotted lines), 2(p^ + MfVlY/^ (solid lines), 175 GeV (dashed fines) and 
(dash-dotted lines) for = A = 175 GeV, as{Mz) = 0.118 and Tt = 1.43. The lower bunch of 
threshold curves (characterized by the peak at around = 348 GeV) is NLO and the upper bunch 
NNLO. Likewise, the lower bunch of high energy curves is 0{as) and the upper bunch 0{al). We 
note that we have not plotted the threshold curves for /i = -y/g^ and the high energy curves not for 
fj, = 25 GeV, which seems to be a rather unnatural choice for each. The formulae for the 0{a1) 
high energy cross section have been taken from Ref. |68]. For convenience we have plotted the curves 
in Fig. 12 in the pole mass scheme. Because the choice of the mass definition does not alter the 
behaviour of the cross section normalization for energies above the peak position, this choice does not 
affect the conclusions drawn below. For the threshold (high energy) cross section, we observe that the 
NNLO {0{ag)) corrections decrease for energies further away from the threshold region. However, 
the 0(0^) corrections to the high energy cross sections are much larger than the NNLO corrections 
to the threshold cross section at the same centre-of-mass energy. For =360-370 GeV the 0{al) 
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Figure 12: The vector-current-induced total cross section in the non-relativistic expansion at NLO 
(lower bunch of threshold curves), NNLO (upper bunch) and in conventional perturbation theory at 
0{as) (lower bunch of high energy curves) and C(ag) (upper bunch). The pole mass scheme has been 
used. The curves have been plotted for OsiM^) = 0.118, Tt = 1.43 GeV, Mt = A = 175 GeV and 
H = 25 GeV (dotted hues), 2{pl^+M^ F^^) V4 (sdid lines), 175 GeV (dashed lines) and (dash-dotted 
lines). The formulae for the 0{a'^) high energy cross section have been taken from Ref. p8| . 



corrections to the high energy cross section are, for equal choices of renormalization scales, between 
10 and 20% compared to only around 5% for the NNLO corrections to the threshold cross section. We 
also see a much weaker renormalization scale dependence of the threshold cross sections. The curves 
show that the resummation of Coulomb singular terms contained in the threshold calculation leads to 
a considerable stabilization of the cross section determined in conventional perturbation theory in ag 
for energies below y/s = 365 GeV. If we believe that conventional perturbation theory is reliable down 
to energies around ^/s = 360 GeV, the results displayed in Fig. 12 indicate that the non-relativistic 
expansion does certainly not break down. However, the curves of Fig. 12 also make it evident that the 
small renormalization scale dependence of the NLO threshold cross section does certainly not reflect 
the true size of the remaining theoretical uncertainties at NLO. We believe that 10% should be a fair 
estimate of the remaining theoretical uncertainties contained in the normalization of the NNLO total 
cross section close to threshold. As far as the top mass determination at a future electron-positron 
linear or muon pair collider is concerned, this rather large normalization uncertainty might in fact lead 
to uncertainties in the determination of the IS" mass that are larger than indicated in the previous 
subsection. This is due to the effects of beamstrahlung and initial state radiation that lead to a 
smearing of the effective centre-of-mass energy of about 1-2 GeV Q . Beamstrahlung and initial state 
radiation render the visible peak in the total cross section either smaller (at the muon pair collider) 
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Figure 13: The total cross section u/of , Eq. (^), is plotted in the 15 scheme at NNLO for Mis = 
A = 175 GeV, Tt = 1.43 GeV and a = 1/125.7. (a) shows the renormalization scale dependence for 
H = 15 (solid line), 30 (dashed line) and 60 GeV (dotted line), and as{Mz) = 0.118; (b) shows the 
dependence on the strong coupling for as{Mz) = 0.113 (solid line), 0.118 (dashed line) and 0.123 
(dotted line), and for /x = 30 GeV. 



or completely invisible (at the linear collider), which makes it possible that the uncertainty in the 
normalization feeds into larger uncertainties in the determination of Mig. It is the task of realistic 
simulation studies to determine how large this effect is for the various collider and detector designs 
and to devise optimized strategies to minimize it. If the effects of beamstrahlung and initial state 
radiation on the top quark mass determination are small, the uncertainty in the normalization will 
mainly affect the measurement of the strong coupling (see Figs. ^ and p!3[b). 

In Figs, m the total cross section cj/of (e+e" ^ 7*, Z* ^ tt), Eq. (H), is plotted at NNLO in 
the IS scheme for Mis = A = 175 GeV, Tt = 1.43 GeV and a = 1/125.7. Figure Ha shows the 
renormalization scale dependence for ^ = 15 (solid line), 30 (dashed line) and 60 GeV (dotted line), 
and as{Mz) = 0.118. Figure |l3| b displays the dependence on the strong coupling for as{Mz) = 0.113 
(solid line), 0.118 (dashed line) and 0.123 (dotted line), and for = 30 GeV. 



9 Summary and Conclusions 

Within the framework of the non-relativistic effective field theories NRQCD and PNRQCD, we have 
calculated the vector-current-induced total cross section of top-antitop pair production in electron- 
positron annihilation close to threshold at NNLO in the non-relativistic expansion. The corresponding 
NNLO QCD relativistic corrections have also been determined for the vector-current-induced top 
three-momentum distribution. In addition, the axial-vector-current-induced total cross section and 
the three-momentum distribution have been calculated to fully account for the Z-boson contributions 
in electron-positron annihilation. For the total cross section and the three-momentum distribution, 
the axial- vector-current-induced contributions are suppressed by with respect to the vector current 
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contributions; they have therefore been determined in leading order in the non-relativistic expansion. 
The size of the axial-vector-current-induced contributions is smaller than the remaining theoretical 
uncertainties in the vector-current-induced cross section (for unpolarized electrons and positrons). In 
contrast with previous literature on the same subject, we have implemented the top quark width by 
including electroweak corrections into the (P)NRQCD matching conditions of the Lagrangian and 
the currents. This allows for a straightforward generalization of the Fadin-Khoze prescription "i? — > 
E + iV^" to implement the top quark width at NNLO in the non-relativistic expansion, where T^/Mf 
is counted as order v'^. We have shown that at NNLO this cannot be achieved by a simple shift of 
the centre-of-mass energy to complex values. Our calculations have been carried out using numerical 
techniques to solve the corresponding integral equations within a cutoff regularization scheme and 
using analytic methods for the matching procedure. We have addressed the question of large NNLO 
corrections to the peak position and the normalization of the total vector-current-induced cross section. 
The position of the peak, which is observable in the total vector-current-induced cross section, can be 
stabilized if the cross section is expressed in terms of the IS mass, instead of the pole mass. The IS" 
mass, Mis, is defined as half the mass of a fictitious ^5*1 toponium ground state for a stable top quark. 
The IS* mass is a short-distance mass and, by construction, reduces to a large extent the dependence 
of the peak position on theoretical parameters such as the renormalization scale of the strong coupling. 
We have also shown that the large NNLO corrections to the normalization of the total cross section, 
of order 20%, are genuine NNLO corrections, which cannot be removed by a redefinition of the top 
quark mass or the strong coupling. The large size of the corrections to the normalization originates 
from the fact that NNLO relativistic corrections from several sources have the same sign. We believe 
that the remaining theoretical uncertainties in the normalization are of order 10%. If the effects 
of beamstrahlung and initial state radiation at the e^e~ linear collider do not lead to a significant 
cross feed of the uncertainties in the normalization into Mis, we expect that an uncertainty in the 
determination of Mis of less than 200 MeV will be possible at the linear collider with an integrated 
luminosity of 50-100 fb^^. In order to determine the MS top quark mass from the IS* mass with the 
same precision, the knowledge of the full three-loop relation between the pole and the MS mass, and 
a small uncertainty in as{Mz) are crucial. 



After completion of this work, we received Refs. |69, 70, O]. In Ref. |p9| the total vector 



current-induced cross section has been calculated analytically, using the MS regularization scheme 
based on the Schrodinger equation (p!9[). The NLO and NNLO corrections have been treated pertur- 
batively, supplemented by a resummation of the energy denominators for the n = 1 and n = 2 states in 
the spectral representation of the Green function. The renormalization scale dependence of the cross 
section line-shape is considerably larger in Ref. |69| than in our work. This might be a consequence 
of the perturbative treatment of the NLO and NNLO corrections. In addition, a next-to-leading log- 
arithmic resummation of logarithms of the ratio Mt/fi in the short-distance coefficient has been 
carried out, taking the MS cutoff scale ^ of order MfV. The effect of this resummation is around 5% 
for the normalization of the total cross section. In our cutoff scheme, where the regularization scale 
is of order Mt, the corresponding logarithm is contained in the non-relativistic current correlators. 
In Ref. [70 1 the vector-current-induced cross total section and the three- momentum distribution have 
been calculated at NNLO, based on the simplified Schrodinger equation (70), which we have discussed 
critically at the end of Sec. ^. For the three- momentum distribution the authors of Ref. |70| have 
included further corrections to account for the difference with the results of the correct Schrodinger 
equation (|55|). In Refs. [^, ^] the "potential-subtracted" mass has been tested in different ways 
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as an alternative mass parameter for the total cross section. As far as the uncertainties in the top 
mass determination at a future linear collider are concerned, Ref. 169(1 arrives at conclusions similar to 



ours. In Ref. |71] the techniques used in Refs. |17, 11| have been employed to calculate the total cross 
section, the angular distribution and the top quark polarization for top quark pair production close to 
threshold in e~^e~ and 77 collisions. The corrections originating from the higher-order contributions 
in the Coulomb potential have been calculated analytically. In Refs. |^, |7l| the top quark width 
has been implemented by the replacement rule ^^E ^ E + iTt" , where E is the centre-of-mass energy 
with respect to two times the top mass. 
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A Calculation of the Short-Distance Coefficient 



In this appendix we present details of the calculation of the short-distance coefficient to order aj, 
assuming that the top quarks are stable (Ft = 0) . We recall that C" is the square of the short-distance 
coefficient cf of the ^Si NRQCD current ip'^irx (see Eq. @); contains those contributions in the 
vector-current-induced top-antitop production diagrams, which come from loop momenta p = (jP,]}) 
with IpI > A for = 2Mt. As explained in Sec. ^ we have to determine C" by employing the 
specific routing convention shown in Fig. |l|. In principle, it would be possible to determine C" by 
calculating the diagrams for the vector-current-induced cross section in full QCD restricting the loop 
momenta such that the spatial components would be larger than A. However, in a cutoff scheme it 
is more economical to first calculate the vector-current-induced cross section in NRQCD up to order 
a1 and NNLO in the velocity expansion and then to adjust the coefficients of C" such that the cross 
section in NRQCD is equal to the cross section in full QCD, likewise calculated to order a1 and NNLO 
in the velocity expansion. 

The expression of the total vector-current-induced cross section in full QCD at order and 
NNLO in the velocity expansion reads (a = Cpas^fJ.)): 
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Figure 14: Graphical representation of the NRQCD vector-current correlators diagrams needed to 
determine the non-relativistic vector-current-induced cross section at the Born level and NNLO in the 
non-relativistic expansion. 




Figure 15: Graphical representation of the NRQCD vector-current correlators diagrams needed to 
determine the non-relativistic vector-current-induced cross section at 0{as) and NNLO in the non- 
relativistic expansion. 
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The Born and 0{as) |7^ contributions are standard. At order the contributions in Eq. ( p9| ) that 
are proportional to Cp, CaCf, CfTui and CfT have been calculated in |7^, |35], |63, 74] and p^, [75[|, 
respectively. (See also Refs. p8|, [7^].) 

To determine the corresponding total vector-current-induced cross section in NRQCD, we 
have to calculate the absorptive part of the correlator diagrams depicted in Figs. 14, |l^ and 0. The 
various symbols are defined in Fig. We emphasize that we neglect multiple insertions of NNLO 
contributions. 

The results for the absorptive parts of the individual diagrams read (a = CfcXs, D{k) = 
Mt/{k^-pl-ie)y. 
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Figure 16: Graphical representation of the NRQCD vector-current correlators diagrams needed to 
determine the non-relativistic vector-current-induced cross section at 0{a^) and NNLO in the non- 
relativistic expansion. 
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Figure 17: Symbols describing the interactions potentials V^°, V^^'^, V^^^°, Vbf and Ina and the 
kinetic energy correction dHj^i^ = (p^ — k'^)/{AMf). V,^'^, V^^° and V^i^nlo pgfgj- i^j-^g Born, one-loop 
and two-loop contributions to the Coulomb potential presented in Eq. (22). 
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where the upper index of the functions /j corresponds to the power of the strong couphng of the 
diagrams and the lower index to the numeration given in Figs. 14, 15 and 16. Combinatorial factors 
are taken into account. We note that the results above have been expanded in pq/Mi^pq/ K <^ 1; no 
condition has been assumed for the ratio A/M^. 

Summing all terms leads to the total vector-current-induced cross section in NRQCD: 
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The short-distance coefficient has to be chosen such that the RHS of Eqs. (|89D and (|107D are equal 
for all terms up to order and NNLO in the non-relativistic expansion. The result reads 
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In Eq. ( |108| ) we have displayed the NLO and NNLO short-distance contributions separately. We 
note that the NLO short-distance contributions in Eq. ( |109| ) differ from the commonly quoted one 
by the term 4^^^. This is a consequence of our cutoff regularization scheme, which excludes 
loop momenta with spatial components larger than A, even if the corresponding integration is UV- 



convergent. We also point out that the NNLO short-distance contribution in Eq. (110) contains the 
term which is of order only. This term is a manifestation of the power-counting breaking 

effects discussed in Sec. ^ The term exists because it subtracts the power-counting breaking terms 
originating from the linear UV-divergent behaviour of the Breit-Fermi potential, Vbf, the kinetic energy 
correction, (pq— (4M^^), and the dimension-5 NRQCD vector-current in the non-relativistic current 
correlators. Thus it is important to consider this term as NNLO. It is a conspicuous fact that the 
NLO short-distance coefficient cJJlo vanishes for the choice A = Mt. We emphasize, however, that this 
cancellation is purely accidental. Nevertheless, comparing the short-distance constant C" calculated in 



our cutoff regularization scheme with the corresponding coefficient obtained in the MS scheme |35, 36] 
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we find that the perturbative corrections are in general smaller in the cutoff scheme. In Table |^ we have 
displayed the NLO and NNLO short-distance corrections for exemplary choices of the renormalization 
scale /X and the cutoff A for as{Mz) = 0.118 and using two-loop running for the strong coupling. 
We are not aware of any principle reason why the short-distance corrections should be, in general, 
better convergent in our cutoff scheme than when using the MS regularization. We finally note that 
the potentially large logarithmic term o^CaCf ln(Mi//x) in C-^, which corresponds to an anomalous 
dimension of the dimension-3 NRQCD vector-current "^^cx, does not exist in C"^, since in our cutoff 
scheme the corresponding logarithmic divergence in the NRQCD diagrams is cut off at the scale 
A ~ Mt rather than ~ Mtv, as in the MS scheme. However, we emphasize that the absence of 
this logarithmic term in C" is traded for the existence of logarithms of the ratio 2pq/A in the non- 
relativistic correlator, which are not present in the MS scheme. The logarithms of A//x and Mt/fi in 
, Eq. (ffTof) , originate from the running of the strong coupling and are not related to an anomalous 
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dimension. 
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0.128 


30 


0.228 
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0.018 


-0.069 


-0.029 


0.025 


60 


0.202 





-0.107 


-0.214 


0.039 


-0.048 


-0.022 


0.037 


175 


0.172 





-0.091 


-0.182 


0.056 


-0.028 


-0.015 


-0.091 



Table 3: The NLO and NNLO contributions to the short-distance coefficient in our cutoff scheme 
and in the MS scheme for various choices of the cutoff A and the renormahzation scale fi. We have 
chosen as{Mz) = 0.118, and two-loop running for the strong coupling has been employed. 
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